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GRATING  INTERFEROMETRIC  SENSORS 

1 

INTRODUCTION 

Conventional  remote  sensing  approaches  provide  Information  on  the 
Intensity  distribution  of  the  object  scene  with  a  resolution  defined  by 
the  real  collection  aperture  of  the  sensor.  In  this  research  project, 
we  explored  the  use  of  coherence  measurements  to  provide  additional 
Information  that  may  help  to  enhance  the  capability  of  a  remote  sensor 
to  discriminate  between  different  targets  or  between  targets  and  the 
background  In  the  visible  wavelengths.  The  enhancement  can  be 
accomplished  by  improving  the  Image  resolution  beyond  that  permitted  by 
the  real  collection  aperture  or  by  providing  information  beyond  the 
Intensity  distribution  of  the  object  scene.  In  particular,  the  use  a 
grating  Interferometer  to  perform  the  coherence  measurements  was 
examined. 

It  Is  useful  to  separate  coherence  measurements  and  their 
applications  into  two  categories:  1)  Measurement  of  the  field  coherence 
at  the  object  plane  where  the  object  scene  is  under  partially  coherent 
illumination.  2)  Measurement  of  the  field  coherence  at  the  detection 
plane  where  the  object  scene  is  self  emissive  or  under  non-coherent 
illumination. 

Measurement  of  the  field  coherence  at  the  object  plane  can  be 
performed  by  mapping  the  object  field  to  the  detection  plane  with  an 
Imaging  system.  We  have  shown  that  for  an  actively  Illuminated  target, 
the  Image  obtained  through  a  scattering  medium  can  be  significantly 
enhanced  by  measuring  the  coherence  of  the  Image  field  with  a  grating 
Interferometer  and  removing  the  contribution  of  tht  incoherent  scattered 
light  to  the  Image  Intensity  [1-3].  Coherence  measurement  at  the  object 
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plane  can  also  be  used  to  obtain  additional  Information  on  a  target 
scene.  When  partially  coherent  field  Is  reflected  by  an  object  surface. 
Its  coherence  properties  may  be  altered  by  differing  amounts,  depending 
on  the  composition  of  the  surface  material  and  texture  (surface  height 
variations).  Thus  by  measuring  the  coherence  of  reflected  field, 
unresolved  surface  height  variations  may  be  Inferred.  Speckles,  which 
are  produced  by  mutual  Interference  of  adjacent  pixels  can  be  used  to 
determine  the  coherence  of  the  object  field.  It  Is  well  known  that  by 
Illuminating  a  surface  with  partially  coherent  light  whose  longitudinal 
coherence  length  Is  of  the  same  magnitude  as  the  surface  height 
variations,  the  texture  of  the  object  surface  can  be  determined  from  the 
contrast  of  the  Image  speckles  [4,5].  Such  a  technique  has  been  applied 
In  metrology  to  study  the  surface  texture  of  machined  parts.  The 
approach  however,  has  not  been  applied  to  study  the  surface  height 
variations  over  a  large  surface  area  for  the  purpose  of  target 
discrimination  In  a  remote  sensing  system. 

If  the  object  field  Is  self-emissive  or  under  non-coherent 
Illumination,  the  field  at  the  object  plane  Is  simply  Incoherent,  no 
useful  Information  can  be  gained  from  measuring  the  coherence  at  the 
object  plane.  However,  there  1$  a  Fourier  transform  relationship 
between  the  Intensity  distribution  of  the  object  scene  and  the  spatial 
coherence  at  a  distance  detection  plane.  This  relationship  Is  utilized 
by  a  new  class  of  passive  sensors  to  achieve  Image  resolution  beyond  the 
diffraction  limits  of  their  Input  apertures.  The  grating  Interferometer 
Is  capable  of  measuring  the  complex  spatial  coherence  of  an  optical 
field  very  efficiently.  We  shall  show  In  this  report,  how  the  grating 
Interferometer  may  be  used  to  detect  the  presence  of  moving  targets  and 
to  perform  passive  synthetic  aperture  Imaging. 

While  we  have  utilized  the  grating  Interferometer  for  object  plane 
coherence  measurements  to  achieve  Improve  seeing  through  scattering 
medium  and  target  discrimination  [1-3],  the  same  approach  can  also  be 
Implemented  with  conventional  Interferometers  albeit  less  easily  In  some 


cases.  The  grating  Interferometer  however,  offers  some  unique 
capabilities  when  used  for  coherence  measurements  at  the  detection  plane 
to  perform  for  example,  synthetic  aperture  Imaging.  For  this  reason,  we 
shall  concentrate.  In  this  report,  on  the  use  of  the  grating 
Interferometer  for  detection  plane  coherence  measurements. 

In  following  sections,  the  van  Cittert-Zernike  theorem  which  forms 
the  basis  for  passive  Interferometric  Imaging  Is  first  reviewed.  The 
measurement  of  the  mutual  coherence  function  with  a  grating 
Interferometer  Is  then  described  and  the  possible  uses  of  the 
Interferometer  for  passive  synthetic  aperture  Imaging  and  as  a  moving 
target  sensor  are  discussed.  An  overview  of  the  grating  interferometric 
sensors  Is  given  In  the  main  body  of  the  report.  The  operation  of  the 
grating  Interferometer  and  Its  uses  In  Incoherent  optical  processing  and 
passive  synthetic  aperture  Imaging  are  described  more  comprehensively  in 
the  articles  [7,8]  attached  as  the  Appendices. 


BASIC  PRINCIPLE 


Astronomers  have  for  some  time  made  use  of  Interferometric 
techniques  for  astronomical  Imaging,  from  Hichelson's  simple 
two-aperture  Interferometer  to  sophisticated  very  long  baseline 
Interferometers  that  span  the  globe.  They  are  all  based  on  the  van 
Cittert-Zernike  theorem  which  relates  the  Intensity  distribution  and  the 
mutual  coherence  function  of  the  fields  In  two  widely  separated  planes. 

Consider  the  system  geometry  Is  Illustrated  In  Figure  1.  An  object 
field  radiates  according  to  Its  blackbody  temperature  and  the  emitted 
radlalton  Is  received  at  (x^.y-j)  and  Assuming  that  the  emitted 

radiation  Is  quasimonochromatic  (or  made  so  by  spectral  filtering),  then 
by  using  the  Huygen-Fresnel  principle,  the  wave  disturbances  at 
and  (012,32)  due  to  the  radiation  from  the  object  field  A  can  be  written 
as  [6]. 
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where  x  Is  the  wavelength  corresponding  to  the  center  frequency  of  the 
filter,  (a|,6|)  and  Eg(a2,B2)  are  the  wave  disturbances  at  (“^,3^) 
and  (a2,32),  and  R2  are  the  distances  between  (01^, 3^)  -  (x-j,  y^ )  and 
(02,^2)  -  (X2,y2)  respectively,  k  Is  the  wave  number  and  the  angular 
frequency  of  the  radiation.  Since  the  radiation  Is  assumed  to  be 
quasimonochromatic,  we  shall  from  here  on,  neglect  the  time  varying 


phase  factor  an  let  E  =  E^e  • 

The  mutual  Intensity  between  Q  and  Q2  Is  given  by 


J(Xpy|;x2.y2)  =  <E(x^  ,y^  )E*(x2,y2)> 


(3) 


where  <->  denotes  ensemble  averaging.  Substituting  Eq.  (1)  and  Eq.  (2) 
into  Eq.  (3),  we  obtain 
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Since  the  emitted  field  is  spatially  incoherent, 
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=  I(*Ri3) 

for 

(5) 

where  I(a,B)  is  the  intensity  distribution  at  the  object  plane.  The 
mutual  intensity  can  therefore  be  simplified  to 


J(Xj^,yj^;x2,y2  )  = 


-ilc(R^-R2) 
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dadiS 


(6) 


The  complex  degree  of  spatial  coherence  (CDSC)  is  related  to  the  mutual 
intensity  by  [6] 
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where 


l(Xj^,yj^)»l(X2,y2) 


=  Ji 


Lt2i£ld.<ia 
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Let  us  assume  that  the  a  -  2  plane  and  the  x-y  plane  are  parallel  to 
each  other,  we  can  then  write 

Rf  =  (x^  -  ot)2  +  (y^  .  s)2  +  r2  (9) 

where  R  is  the  separation  between  thetarget  plane  and  the  detection 
plane. 


R-j  can  be  approximated  by 
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and  similarly 


Substituting  Eq.  2-14  Into  Eq.  2-7  we  obtain 
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where  u  =  (x^  -  X2)/RX,  v  =  (y^  -  y2)/RX  and 

2  2  2  2 

(xj+yi)-(x2‘^y2^ 

“  2RA 

The  CDSC,  ’•*  equal  to  the  Fourier  transform  of  the  object 

Intensity  distribution  multiplied  by  a  quadratic  phase.  If  the  far 
field  condition  Is  satisfied,  then  =  0  and  e^'^  =  1.  The  CDSC  becomes 
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where  K  Is  a  complex  constant. 

We  see  that  the  CDSC  Is  related  to  the  Intensity  distribution  of 
the  object  field  by  a  Fourier  transform  relationship.  We  can  therefore 
reconstruct  the  object  Image  by  measuring  the  two-dimensional  coherence 
function  and  performing  an  Inverse  Fourier  transform. 
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The  COSC  can  be  obtained  with  Young's  simple  two-pinhole 
arrangement  or  equivalently,  Michelson's  stellar  Interferometer,  by 
measuring  the  complex  visibllty  of  the  Interference  fringes.  The 
complex  fringe  visibllty  function  Is  defined  as 


^l(x  ,y  ) 


Except  for  very  near  field  conditions,  we  can  as  before  make  the 
approximation,  I(x^,y2)  »  I(x2,y2)  and  we  have 


V(u,v)«r  y(u,v) 


To  measure  the  COSC  via  the  visibility  function  Is  not  very  efficient 
since  the  Interference  pattern  must  be  scanned  out  with  a  small  detector 
to  determine  the  phase  and  modulation  depth  of  the  sinusoidal  fringe 
pattern.  An  alternate  approach  Is  to  use  a  beam  splitter  to  recombine 
the  fields  gathered  at  (x^,y|)  and  propagate 
coincidentally.  The  entire  output  can  then  be  focussed  on  and  detected 
by  a  photo  detector.  It  can  be  shown  that  the  output  corresponds  to  the 
Inphase  portion  of  the  Fourier  transform  sampled  at  spatial  frequency 
(u,v)  (see  Appendix  8).  The  quadrature  part  can  be  obtained  by  shifting 
the  position  of  the  beam  splitter  to  generate  a  -n/Z  phase  difference 
between  the  recombined  fields.  By  changing  the  positions  of  (x^,y^)  and 
(X2>y2)>  different  spatial  frequencies  of  the  object  Is  measured.  A 
large  aperture  can  then  be  synthesized  In  time. 

We  can  see  from  Eq.  15  that  the  Fourier  transform  kernel  Is  scaled 
by  wavelength.  This  fact  Is  used  In  synthetic  aperture  Imaging  systems 
to  generate  the  radial  fill  of  the  synthetic  aperture.  This  wavelength 
dependency  however,  also  limits  the  Instantaneous  spectral  bandwidth 


•  »  *  - 


•  V  A 


that  can  be  used  with  either  the  active  or  the  passive  systems.  Opening 
the  bandwidth  too  wide  will  produce  a  smearing  effect  In  the  phase 
history  signal  and  reduce  the  system  signal-to-nolse  ratio.  This 
restriction  does  not  pose  any  problem  for  active  systems  but  It  limits 
the  performance  of  passive  Interferometric  Imaging  systems.  Natural 
emissions  and  Illuminating  sources  are  spectrally  very  wide  band  and  the 
SNR  of  the  Images  reconstructed  with  a  passive  Interferometric  Imaging 
system  Is  directly  proportional  to  the  Instantaneous  spectral  bandwidth. 
System  sensitivity  can  therefore  be  Improved  If  a  means  can  be  found  to 
Increase  the  Instantaneous  spectral  bandwidth  without  suffering  the 
signal  loss  due  to  the  wavelength  dependency  of  Fourier  transformation. 

In  the  next  section  the  use  of  an  achromatic  grating  Interferometer 
for  passive  synthetic  aperture  Imaging  of  rotating  objects  In  the 
visible  region  Is  described.  The  achromatl cl ty  of  the  grating 
Interferometer  permits  the  Instantaneous  spectral  bandwidth  to  be 
widened  up  substantially  while  maintaining  the  capability  to  use 
wavelength  diversity  to  achieve  aperture  fill.  And  as  we  shall  show* 
the  achromatic  grating  Interferometer  also  offers  other  operational 
advantages  over  conventional  Interferometric  systems  when  used  for 
passive  synthetic  aperture  Imaging. 
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ACHROMATIC  GRATING  INTERFEROMETER 

The  achromatic  grating  Interferometer  has  been  studied  extensively 
and  applied  to  many  applications  [1,8-10].  Tal  and  Aleksoff  In 
particular,  has  utilized  It  as  an  Incoherent  optical  processor  for  real 
time  Fourier  transformation  [10].  The  transform  operation  Is  limited  to 
a  single  dimension  and  It  could  create  difficulties  in  some 
applications.  This  limitation  however,  does  not  represent  any  drawback 
for  synthetic  aperture  Imaging  since  under  most  operating  conditions,  a 
long  baseline  Is  available  only  along  a  single  dimension  (e.g.,  along 
the  fuselage  of  an  airplane). 

Since  the  Interferometer  performs  the  correlation  only  along  a 
single  dimension  at  any  one  time,  a  two  dimensional  model,  comprising  of 
the  direction  of  propagation  and  the  direction  of  lateral  separation  of 
the  fields  to  be  correlated,  will  be  used  to  describe  the  operation  of 
the  grating  Interferometer.  In  Figure  2,  we  show  that  geometry  of  a 
grating  Interferometer  that  can  be  used  for  passive  synthetic  aperture 
Imaging.  It  Is  composed  of  three  gratings  of  the  same  spatial  frequency 
F.  Light  waves  Impinging  on  gratings  G^  and  6^  are  diffracted  and 
recombined  at  G^  which  Is  placed  a  distance  d  from  the  plane  formed  by 
and  G2.  The  gratings  G^  and  also  represents  the  Input  apertures 
and  they  have  a  width  of  A  and  are  separated  by  a  distance  S.  After 
being  diffracted  by  G^,  the  recombined  fields  propagate  coincidentally 
and  are  focussed  down  onto  a  photo  detector  of  width  W^.  Let  us  assume 
as  before  that  the  object  plane  Is  sufficiently  far  away  that  the  light 
emanated  by  a  point  source  In  the  object  scene  may  be  considered  to  be  a 
plane  wave  as  It  impinges  on  the  Input  apertures  of  the  grating 
Interferometer.  Thus,  for  a  single  point  source  of  Intensity  OCa^j) 
located  at  a  distance  from  the  optical  axis  as  shown  In  Figure  3,  the 
optical  field  just  before  It  Is  diffracted  by  the  grating  G3  Is  equal  to 


E(«o)  *  ^  ~  -  ikd  ^  1  -  x'^Ip  -  F*")] 

+  exp  [i2ir(p  +  F)x  -  ikR'/~l  -  x^p^  -  ikd  /  1  -  x^(p  +  F)^]  |  (17) 

where  F  is  the  spatial  frequency  of  the  diffraction  gratings,  p  = 
sin[atn(oi^/R)]/X  and  k=2T:/X.  The  corresponding  intensity  distribution  can 

be  written  as 

Ig(a  )  »  C0{a^)  I  1  +  cos  [4,tFx  -  -  x^lp  +  F)^  -  /l  -  x^Cp  -  F)^  ) 

(18) 

where  C  is  an  appropriate  constant.  Using  a  grating  63  with  spatial 
frequency  F  to  recombine  the  beams  diffracted  by  61  and  62,  the  field 
intensity  is  demodulated  back  to  base  band  in  terms  of  spatial  frequency 
and  the  intensity  distribution  right  after  being  diffracted  by  63 
becomes 

Io(c«o)  *  CO(aQ)  j  1  -  cos[i^  (/l  -  x^(p  +  F)^  -  /  1  -  x2(p  -  F)^  )]| 

(19) 

If  we  expand  the  square-root  terms  into  binomial  series  and  keep  only 
the  first  order  terms,  we  obtain 

2Fda  1 

iQ(ao)  =  CO(ag)  [  1  -  COS  2n(— ^)J  (20) 


We  see  that  the  output  intensity,  to  the  first  order,  is  not  a  function 
of  wavelength.  The  object  scene  within  the  FOV  of  the  interferometer 
can  be  considered  as  a  superposition  of  incoherent  point  sources.  The 
output  intensity  due  to  the  entire  scene  can  then  be  expressed  as 


where 


da  +  C  y*  0(a) 
Z 
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u 


2Fd 

•IT 
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which  gives  a  cosine  transform  relationship  together  with  a  bias  term. 
If  we  translate  grating  G3  laterally  by  1/8  of  a  grating  period,  the 
phases  of  the  two  fields  being  recombined  are  shifted  achromatically  by 
^^/4.  It  can  be  shown  that  the  output  Intensity  becomes 


H  C  /o(a) 
Z 


da  +  C 


/o( 


a)  s1n  (2vUa)  da 


(22 


giving  the  sine  transform  of  the  Intensity  distribution  of  the 
Incoherent  source. 

For  a  given  value  of  F,  d  and  R,  the  grating  Interferometer 
measures  the  source  spectrum  at  a  spatial  frequency  of  u  =  Fd/R.  Unlike 
the  conventional  system  described  earlier,  the  transform  kernel  Is  not 
wavelength  dependent. 

A  more  detailed  discussion  of  the  grating  Interferometer  and  Its 
achromaticity  are  given  In  Appendices  A  and  B. 


4 

PASSIVE  SYNTHETIC  APERTURE  IMAGING 

For  a  given  Interferometer  geometry,  the  correlation  output 
provides  the  value  of  the  object  spectrum  at  a  single  spatial  frequency. 
The  system  geometry  must  be  changed  to  generate  transform  data  at  other 
spatial  frequencies  In  order  to  synthesize  a  two-dimensional  aperture. 
First,  If  there  Is  relative  rotational  motion  about  the  optical  axis 
between  the  object  scene  and  the  Interferometer,  transform  data  can  be 
gathered  along  a  circular  arc  In  the  transform  plane.  Second,  If  either 
the  baseline  or  the  detection  wavelength  Is  varied,  transform  data  long 
a  radial  line  In  the  transform  plane  Is  obtained.  Combining  these  two 
effects,  a  two-dimensional  aperture  can  be  synthesized.  The  means  to 
achieve  such  an  aperture  fill  will  be  discussed  below. 

Let  us  assume  for  now  that  the  Interferometer  Is  looking  straight 
down  at  the  object  scene  located  at  the  optical  axis.  The 
Interferometer  output  provides  the  transform  data  at  a  single  point  In 
the  Fourier  transform  plane.  As  the  object  (or  equivalently,  the 
Interferometer)  rotates  about  the  optical  axis,  data  Is  gathered  along  a 
circular  arc. 

Radial  fill  of  the  synthetic  aperture  can  be  achieved  with  the  use 
of  baseline  diversity  or  wavelength  diversity.  Since  the  spatial 
frequency  measured  Is  determined  by  S/RX,  It  can  be  changed  by  varying 
S,  the  separation  of  the  detection  aperture  (baseline)  or  X,  the 
detection  wavelength.  Baseline  diversity  can  be  accomplished  with  the 
grating  Interferometer  by  moving  the  three  gratings  closer  or  apart. 
However,  It  may  not  the  desired  mode  of  operation  for  synthetic  aperture 
Imaging  since  It  requires  the  use  of  large  Input  windows  to  accomodate 
the  change  In  baseline.  It  Is  often  more  preferrable  to  keep  the  Input 
aperture  fixed  and  achieve  the  radial  aperture  via  the  use  of  wavelength 
diversity.  On  the  surface.  It  may  seem  that  the  use  of  wavelength 
diversity  with  an  achromatic  system  Is  contradictory.  We  note  however, 
that  the  Interferometer  has  a  finite  Instantaneous  bandwidth  and  the 


whole  band  can  be  sllded  back  and  forth  to  obtain  the  wavelength 
diversity  for  radial  aperture  fill.  This  Is  accomplished  by  changing  d. 
the  separation  between  gratings  61-G2  and  G3.  Recall  from  Eq.  24  that 
the  spatial  frequency  measured  by  the  Interferometer  Is  equal  to  2Fd/R. 
Changing  d  thus  varies  the  spatial  frequency  measured.  It  can  be  seen 
from  Figure  3  that  changing  d  Is  also  equivalent  to  changing  the  center 
wavelength  of  measurement  X. 

The  two  means  of  aperture  (circular  and  radial)  fill  can  be  used 
together  to  generate  a  two-dimensional  aperture.  Unlike  tomographic 
Imaging,  synthetic  aperture  Imaging  systems  generally  gather  transform 
data  over  an  angular  range  much  less  than  90  degrees  and  the 
two-dimensional  aperture  synthesized  Is  typically  a  segment  of  the  an 
annulus  as  shown  In  Figure  4.  The  Image  (half-power)  resolution  In  the 
ground  plane  Is  approximately  equal  to 


.  Rx  Rx^ 

Rx  Rx 

^  ■  2  s  si7T4" 


(23) 


where  A9  Is  the  angle  of  rotation.  The  resolution  In  the  a  direction  Is 
2 

proportional  to  X  Instead  of  X .  It  Is  because  the  size  of  the 
synthetic  aperture  Is  proportional  to  AX/X,  the  aperture  Is  larger  for 
smaller  X  if  AX  is  fixed.  For  example,  a  real  aperture  Imaging  system 
operating  at  0.5  m  with  a  10  cm  aperture  can  achieve  a  ground  plane 
resolution  of  5m  x  5m  at  a  range  of  1000  km.  A  passive  synthetic 
aperture  Imaging  system  with  the  same  real  aperture  size  and  an  aperture 
separation  of  2m  on  the  other  hand,  can  achieve  a  ground  plane 
resolution  of  0.63m  x  0.58m  at  the  same  distance  with  X  »  0.5  urn,  AX  » 
0.2  urn,  and  A9  »  25®. 


We  have  shown  (see  Appendix  6)  that  a  substantially  wider 
Instantaneous  spectral  bandwidth  can  be  used  with  the  grating 
Interferometer  than  with  a  conventional  Interferometer,  thereby  giving 
the  grating  Interferometer  better  sensitivity.  There  are  also 
operational  advantages  In  using  the  grating  Interferometer.  In  a 
conventional  system,  the  measurement  wavelengths  and  thus  the  spatial 
frequencies,  are  selected  with  a  bank  of  narrow  band  spectral  filters. 
With  the  grating  Interferometer  on  the  other  hand.  It  Is  determined  by 
the  position  of  the  grating  which  can  be  varied  continuously. 
Therefore,  continuous  radial  fill  Is  difficult  to  achieve  with 
conventional  Interferometers  (without  changing  the  baseline)  but  simple 
to  Implement  with  a  grating  Interferometer.  Due  to  the  sampling 
requirements,  measurement  at  many  wavelengths  Is  often  required  to 
generate  an  adequate  aperture  fill.  The  number  of  parallel  Chanels 
required  could  be  fairly  large  (>30).  The  grating  Interferometer 
provides  the  option  of  using  a  single  or  a  small  number  of  parallel 
channels.  The  aperture  fill  Is  obtained  by  simply  dithering  a  single  or 
a  small  set  of  gratings. 

The  grating  Interferometer  Is  also  less  sensitive  to  system 
vibration  and  accurate  phase  shifts  can  be  more  easily  achieved.  It  Is 
because  the  phases  of  a  diffracted  waves  are  determined  mainly  by  the 
lateral  position  of  the  grating  fringes.  Moving  a  grating  laterally  by 
one  grating  period  Instead  of  one  wavelength  produces  one  full  wave  of 
phase  shift.  It  Is  therefore  many  times  less  sensitive  to  component 
motion  than  conventional  interferometers  that  employ  mirrors  and  beam 
splitters. 

A  more  comprehensive  discussion  of  the  use  of  the  grating 
Interferometer  for  passive  synthetic  aperture  Imaging  Is  given  In 
Appendix  B. 
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5 

MOVING  TARGET  SENSOR 

In  some  applications,  the  ability  of  an  Interferometric  sensor  to 
measure  a  target  spectrum  at  high  spatial  frequencies  can  be  used 
without  the  need  to  build  a  full  two-dimensional  aperture.  One  such 
application  Is  the  detection  of  a  moving  unresolvable  point  target. 

Moving  target  detection  can  be  performed  by  various  means  such  as 
the  subtraction  of  time  delayed  Images.  These  techniques  however.  Is 
effective  only  If  the  target  Is  of  a  resolvable  size.  The  energy  of  an 
unresolved  point  object  Is  spreaded  over  the  resolution  cell  of  the 
Image  and  the  resulting  signal-to-nolse  ratio  may  not  be  adequate  for 
unambiguous  target  detection.  The  SNR  can  be  Improved  by  high  pass  or 
bandpass  filtering  since  the  target-to-background  energy  ratio  1s  larger 
at  high  spatial  frequency  than  near  dc  for  a  point  target.  An 
Interferometer  with  a  long  base  line  can  be  used  to  detect  the  presence 
of  a  moving  point  target  that  Is  too  small  to  be  resolved  by  the  real 
Imaging  aperture.  Once  again,  the  need  to  operate  with  a  narrow 
spectral  bandwidth  restricts  the  efficiency  of  conventional  systems  when 
used  1n  a  passive  mode. 

Consider  the  grating  Interferometer  depicted  earlier  in  Figure  2. 
The  output  Intensity  as  a  function  of  target  position  was  found  to  be 
approximately  equal  to 

Ig(u)  =  C0(a^)[l  +  cosCZiTua^j)]  (24) 

where  C  Is  the  constant,  0(01^)  Is  the  Intensity  of  the  point  object  and 
u  »  2Fd/R.  for  a  point  target  moving  linearly  along  the  a  direction, 
the  output  Intensity  becomes 

lQ(u,t)  =  CO(a^)  [1  +  cos(2TTu[ajj  +  Tt)]]  (25) 
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where  T  is  the  translation  rate  and  we  have  assumed  that  the  intensity 
of  the  point  target  does  not  change  with  position  or  time.  A  linearly 
moving  target  thus  produces  a  sinusoidal  output  intensity  which  can  be 
filtered  to  enhance  the  SNR.  The  achromatici  ty  of  the  grating 
interferometer  once  again  leads  to  better  system  sensitivity  by  allowing 
the  use  of  a  wider  instantaneous  spectral  bandwidth. 
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SUMMARY 

-v. 

We  have  examined  the  uses  of  a  grating  Interferometer  for  passive 
synthetic  aperture  Imaging  to  achieve  Imaging  resolution  beyond  that 
permitted  by  the  collection  aperture  and  as  a  moving  target  sensor  to 
detect  the  presence  of  moving  targets  that  are  too  small  to  be  resolved. 
We  have  shown  that  the  achromaticity  of  the  grating  Interferometer 
allows  the  use  of  a  wider  Instantaneous  spectral  bandwidth  which  In  turn 
provides  a  better  system  sensitivity  than  conventional  systems.  In 
addition,  the  grating  interferometer  offers  a  simple  means  to  achieve 
two-dimensional  aperture  fill  and  It  Is  more  tolerant  to  system 

vibration. 
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Grating-based  interferometric  processor  for  reai-time 
opticai  Fourier  transformation 


Anthony  M.  Tai  and  Carl  C.  Aleksoff 


A  processing  approach  is  introduced  that  is  capable  of  performing  1-D  real-time  Fourier  transformations 
on  the  intensity  distribution  of  an  incoherent  optical  input.  The  processing  approach  is  based  on  grating 
interferometers,  and  the  resulting  processors  are  simple  in  structure  and  easily  implemented.  Possible  pro¬ 
cessor  configurations  together  with  experimental  resulta  demonstrating  the  operation  of  the  system  are  pre¬ 
sented.  Analyses  are  given  comparing  the  grating  interferometric  processor  to  the  Michelson  stellar  inter¬ 
ferometer  and  the  classical  coherent  optical  processor. 


I.  Iniroductiofi 

The  Fourier  transformation  is  a  basic  operation  in 
many  signal  processing  applications,  and  it  has  been 
implemented  by  various  optical  techniques.  Coherent 
optical  processors^  can  perform  Fourier  transformation 
in  real  time,  but  they  have  several  familiar  limitations. 
Most  coherent  optical  processing  schemes  require  the 
use  of  an  incoherent-to-coherent  converter  which  re¬ 
duces  the  system  speed,  dynamic  range,  and  linearity. 
An  achromatic  Fourier  transform  system^  s  can  be  uti¬ 
lized  to  relax  the  temporal  coherence  requirement  and 
improve  the  system  SNR  by  reducing  coherence  noise. 
The  spatial  coherence  requirement,  however,  is  un¬ 
changed,  and  an  incoherent-to-incoherent  converter  is 
still  needed. 

One-dimensional  processing  using  bulk  acoustical 
modulators  or  SAW  devices,^  ®  on  the  other  hand,  are 
often  too  limited  in  their  space-bandwidth  product  by 
physical  constraints  such  as  transducer  frequency  re¬ 
sponse,  acoustical  velocity,  and  acoustical  attenuation. 
Moreover,  processing  with  an  acoustical  modulator  is 
limited  to  serial  inputs.  To  be  able  to  process  parallel 
spatial  inputs  such  as  images,  a  spatial-to-temporal 
conversion  is  needed. 

Incoherent  optical  processing  techniques  such  as 
OTF  synthesis®  '^  permit  the  use  of  incoherent  optical 
inputs,  but  their  applications  are  limited  to  a  rather 
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restrictive  class  of  operation.  For  example,  complex 
Fourier  transformation  has  not  been  demonstrated. 
Cosinusoidal  transformation  was  demonstrated  by 
Mertz  and  Rogers®  using  a  shadow  casting  technique. 
Such  a  technique  operates  in  the  geometric  optics  re¬ 
gime,  and  the  achievable  space-bandwidth  is  conse¬ 
quently  snudl.  The  optical  vector-matrix  multiplier®  *® 
can  be  configured  to  perform  Fourier  transformation; 
unlike  other  more  conventional  optical  processors,  the 
vector-matrix  multiplier  operates  on  discrete  sampled 
signals.  The  difficulties  in  constructing  and  aligning 
a  large  2-D  mask  and  the  limits  on  the  densities  of  these 
masks  restrict  the  space-bandwidth  product  that  can 
be  achieved  in  practice.  The  relative  complexity  of  the 
system  also  diminished  its  attractiveness. 

In  this  paper,  we  introduce  an  optical  processor  which 
belongs  to  a  class  of  processor  which  we  call  interfero¬ 
metric  processors.  This  processor  operates  on  inco¬ 
herent  optical  inputs,  but  its  performance  is  closely 
related  to  coherent  optical  processing  systems.  An 
achromatic  grating  interferometer  is  used  in  our  im¬ 
plementation  of  the  processing  approach.  It  is  capable 
of  producing  and  displaying  in  real  time  the  real  and 
imaginary  parts  of  the  Fourier  transform  of  a  1-D  in¬ 
coherent  optical  input.  Preliminary  experimental  re¬ 
sults  are  presented  to  verify  the  theory  and  to  demon¬ 
strate  the  processing  approach. 

II.  Comptox  D«gr*«  of  Spatial  Coharanca 

Consider  the  system  geometry  illustrated  in  Fig.  1. 
The  radiation  emitted  by  a  space-limited  object  field 
at  the  (n,d)  plane  is  detected  at  far-field  locations 
(xi,>'i)  and  (x>,>'2).  If  the  radiation  is  quasi-mono- 
chromatic  and  spatially  incoherent,  it  can  be  shown  that 
the  complex  degree  of  spatial  coherence  (CDSC)  is  equal 
to" 


Pig.  1.  System  geometry  relating  the  intensity  distribution  of  an 
incoherent  object  Held  to  the  complei  degree  of  spatial  coherence  at 
a  detection  plane  in  the  far  field. 


'(u.i  l  =  K  /la.dtexp  luo  cd)  dadd.  (11 


where 


u  =  (jti  -  xiMR.  i'  »  (yi  -  y‘2)/R. 

and  K  is  a  constant. 

We  see  that  there  is  a  Fourier  transform  relationship 
between  the  input  intensity  distribution  and  the  CDSC 
at  far  field.  This  relationship  is  sometimes  known  as 
the  van  Cittert-Zernike  theorem.  Thus,  to  obtain  the 
Fourier  transform  of  the  intensity  distribution  of  an 
incoherent  input,  one  merely  has  to  devise  a  means  to 
measure  the  CDSC  at  far  field.  This  can  be  achieved 
interferometrically  by  measuring  the  complex  visibility 
of  the  interference  fringes  form^  between  the  fields  at 
(jci.v'i)  and  (xj.y.').  The  fringe  visibility  function  is 
equal  to 


V'1(X|.>  il.lr.>,vj)| 


-V  /(.f|,,V|)  V  /(x.>..Vj> 
/(xi.yil  +  /(x.>,_v...) 


•  (’((.ri.yi),(Xi.y2l|.  (2) 


For  a  homogeneous  incoherent  input.  /(xi,yi)  = 
/ (.V  j.v  j),  and  the  CDSC  is  space  invariant  in  the  far  field. 
Thus  we  may  write 


V'lZ/.Cf  *  Cll/.C) 
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With  the  Michelson  stellar  interferometer  and  other 
similar  interferometric  imaging  systems,  the  spatial 


frequencies  are  sampled  by  physically  varying  the 
separation  of  the  detecting  apertures.  Such  a  sequen  - 
tial  measurement  scheme  is  too  restrictive  for  use  in  a 
real-time  processing  system.  The  requirements  that 
the  input  be  quasi-monochromatic  also  makes  the  sys¬ 
tem  very  inefficient  for  inputs  that  are  naturally  poly¬ 
chromatic. 

Roddier  et  al.  ‘  -  introduced  a  rotation  shearing  in¬ 
terferometer  for  astronomical  imaging.  Basically  the 
input  wave  front  is  separated  into  two  parts.  One  part 
is  rotated  180°  with  a  system  of  roof  prisms  or  dove 
prisms  and  then  recombined.  At  the  center  of  the  field 
which  is  the  center  of  rotation  there  is  no  shear.  Mov¬ 
ing  away  from  the  center,  increasing  amounts  of  shear 
are  introduced.  Thus,  unlike  the  Michelson  stellar 
interferometer,  all  the  spatial  frequencies  can  be  dis¬ 
played  at  once.  This  technique  is  attractive  for  many 
applications,  and  it  was  suggested  as  an  alternative  to 
stellar  speckle  interferometry.  However,  similar  to  the 
Michelson  stellar  interferometer,  it  requires  quasi¬ 
monochromaticity  and  is,  therefore,  rather  inefficient 
in  the  processing  of  inputs  that  are  spectrally  wide  band. 
George  and  Wang‘S  combined  the  rotation  shearing 
interferometer  with  an  achromatic  transform  optical 
system  and  demonstrated  polychromatic  cosine  trans¬ 
form.  Their  system,  however,  requires  a  system  of 
prisms  to  perform  the  rotation  and  a  multielement  op¬ 
tical  system  for  the  achromatic  transformation. 
Moreover,  switching  from  cosine  to  sine  transformation 
cannot  be  easily  achieved  due  to  the  difficulty  in  real¬ 
izing  achromatic  ir/2  phase  shifts  using  reflective  and 
refractive  optics. 

In  this  paper,  we  introduce  a  grating-based  achro¬ 
matic  optical  processor.  It  is  capable  of  producing,  in 
real  time,  cosine  and  sine  transformations  of  inputs  that 
are  spatially  and  temporally  incoherent.  One  major 
advantage  offered  by  grating-based  systems  is  sim¬ 
plicity.  The  grating  interferometric  processor  tends  to 
require  much  less  hardware  than  conventional  systems 
making  it  easier  to  fabricate  and  set  up.  The  simplicity 
of  the  grating  system  also  makes  the  processor  more 
compact  and  vibration  resistant. 

III.  Achromatic  Grating  Intarforometar 

The  achromatic  grating  interferometer  has  been 
studied  extensively  and  applied  to  many  applications. 
Chang  and  Leith*^*‘®  applied  the  interferometer  to  in¬ 
terferometric  imaging  and  nondestructive  testing. 
Leith  and  Roth' "  studied  the  noise  performance  of  the 
interferometer  and  the  synthesis  of  the  convolution 
integral  for  linear  filtering.  Collins'®  demonstrated  the 
construction  of  holograms  and  matched  filters  with 
incoherent  light.  Leith  and  Swanson^"  utilized  the 
interferometer  for  the  fabrication  of  noise-free  dif¬ 
fractive  optical  elements.  Tai  and  AleksofP*  employed 
the  grating  interferometer  arrangement  in  an  imaging 
system  for  improved  imaging  through  scattering 
media. 

Leith  and  Chang,"’  in  particular,  have  recognized  the 
ability  of  the  grating  interferometer  to  perform  Fourier 
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Fig.  2.  Single-grating  achromatic  interferometer:  (a)  basic  geom¬ 
etry;  (b)  using  a  beam  splitter  and  a  telescopic  imaging  system  to 
decode  the  processor  output;  (c)  imaging  the  fringe  volume  to  a  plane 
away  from  the  grating  to  access  the  zero  spatial  frequency. 


transformations.  Their  interest,  however,  was  pri¬ 
marily  in  imaging  where  the  image  of  an  object  was  re¬ 
constructed  by  performing  an  optical  transform  with 
a  coherent  optical  processor.  To  be  able  to  utilize  the 
interferometer  as  a  processor,  an  appropriate  decoding 
technique  must  be  developed  to  produce  a  complex 
output  that  can  be  efflciently  read  out  or  interfaced  to 
a  computer. 

The  simplest  achromatic  grating  interferometer  that 
can  be  realized  is  composed  of  a  single  grating.  Con¬ 
sider  the  optical  arrangement  depicted  in  Fig.  2(a).  For 
a  single  point  source  located  at  a  >  ao,  the  field  ampli¬ 
tude  in  the  volume  where  the  diffracted  orders 
overlap  is  equal  to 

Eioupij)  “  Cy/S(ao) 

•  |exp|i2r(Pa  -  F)x  +  i*zv/l  -  XHPo  -  F)’) 

+  expli2T(Po  +  F)i  +  ikzy/ 1  -  X*(Po  +  P)’|l.  (4) 

where  k  *  2ir/\,S(ao)  is  the  intensity  of  the  point 
source,  C  is  a  constant,  and  Pq  »  sin(tan“*ao//)/A,  f  is 
the  focal  length  of  the  lens,  F  is  the  spatial  frequency  of 
the  grating,  and  X  is  the  wavelength  of  the  source,  llie 
corresponding  light  intensity  distribution  in  the  fringe 
volume  can  be  written  as 

/(ao;x^}  •  |C|*S(ao)  |l  +  co»  j4irP*  + 

•  Iv  1  -  X*(Po  +  /')’-  V 1  -  X’(Po  -  F)''‘||  •  (51 


Using  a  first-order  approximation,  the  intensity  dis¬ 
tribution  can  be  expressed  as 


■  |C|25(ao)  |l  +  co« 


4irP  X  -  •  (6) 

/  / 


Now,  if  instead  of  having  a  single  point  source,  we  have 


a  spatially  incoherent  line  source  with  intensity  distri¬ 
bution  /(a)  extending  from  a  =  0  to  a  »  /,  the  output 
intensity  distribution  within  a  constant  factor  is 


/(X^>' 


1  ra\ 

1  f(a) 

1  +  cos4irF  X  — - 

Jo 

\  f 

da. 


(7) 


The  intensity  along  the  optical  axis  at  x  =  0  is 


/(old  +  co»ua)da 


=  r  f(a)da  +  r  /(a)  coiuada.  (81 

Jo  Jo 

where  u> »  AirFz/f. 

The  first  term  is  simply  a  constant  bias,  and  the  sec¬ 
ond  term  we  recognize  as  the  cosine  transform  integral. 
Thus  by  reading  out  the  intensity  distribution  along  the 
optical  axis,  we  obtain  the  bipolar  cosine  transform  of 
the  1-D  input  intensity  distribution  plus  a  constant  bias. 
We  emphasize  that  the  expression  in  Eq.  (8)  is  not  a 
function  of  wavelength.  The  transformation  is  per¬ 
formed  achromatically. 

The  single-grating  interferometer  provides  an  ex¬ 
tremely  simple  means  of  achieving  optical  transfor¬ 
mation  achromatically.  However,  the  information  is 
displayed  as  an  intensity  distribution  along  the  optical 
axis,  which  is  rather  inaccessible.  In  particular,  the  zero 
spatial  frequency  is  located  at  z  =  0,  right  at  the  surface 
of  the  grating.  One  way  to  read  out  the  intensity  dis¬ 
tribution  along  the  z  axis  is  to  place  a  very  thin  beam 
splitter  along  the  z-y  plane  at  x  >  0  as  shown  in  Fig. 
2(b).  The  intensity  distribution  can  then  be  imaged  to 
a  convenient  output  plane  to  be  measured.  Unfortu¬ 
nately,  it  is  not  straightforward  to  place  such  a  beam 
splitter  right  up  against  the  grating.  It  nevertheless  can 
be  done  using  special  optics.  One  may  image  the  fringe 
volume  directly  with  a  telescopic  imaging  system  as 
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Fig. Multiple-grating  achromatic  interferometer. 


shown  in  Fig.  2(c)  to  make  the  fringe  volume  more  ac¬ 
cessible.  However,  the  diffraction  angles  would  be  se¬ 
verely  limited  by  the  numerical  apertures  of  the  lenses. 
A  better  approach  is  to  utilize  a  three-grating  arrange¬ 
ment  similar  to  that  of  Leith  and  Chang'^  as  shown  in 
Fig.  3.  The  first  grating  G  i  has  a  spatial  frequency  of 
F,  the  second  and  third  gratings  are  coplanar,  and  both 
have  a  spatial  frequency  of  2F.  The  input  light  field  is 
split  into  two  first  orders  and  the  first  grating  and  then 
recombined  by  the  second  and  third  gratings,  G2  and 
G3.  A  new  fringe  volume  is  formed  about  z  =  2D.  ac¬ 
cessibly  away  from  any  grating.  Moreover,  spurious 
diffraction  orders  can  now  be  easily  removed  to  elimi¬ 
nate  the  spurious  fringes  and  added  bias  that  may 
otherwise  be  introduced.  With  a  single  point  source 
located  at «  *  ooi  the  light  field  at  the  new  fringe  vol¬ 
ume  is  equal  to 

£(aoU.\  I  «  Cv  5(«o)  •  lexp((2ir(Ai  -  fix 
iitD\  1  —  A*(Po  -f  F)^ 

+  ik[z  +  Olv  1  —  X-(Po  “  f ) *1 
+  exp|i2ir(Po  fix  +  ikD\  1  -  X-(Po  -  f  1^ 

■t-iktz  -t-DIv  1  -  I9I 

The  corresponding  light  intensity  distribution  can  be 
written  as 


/(X.X)  »  Id^SlaKl  -*•  cos|4irfx 

+  *x|v  1  -  \-(Po  +  F)^  -  V  1  -  y^Pa  -  fdi).  (10) 

which  can  be  approximated  by 


I{X.Z)  ■  |Cr-cS(a)  |l  +  COS 


4irF(x-£^]| 


(III 


For  a  1-D  signal  along  the  «  axis  ranging  from  0  to  (. 
the  output  intensity  distribution  becomes 


/Ix.xl  =  ^  /(iiljl  +ci)sj4irF|x  — ^jjjdci  Il2l 

At  X  =0.  once  again  we  have  the  cosine  transform  re¬ 
lationship 


To  perform  a  Fourier  transformation,  both  real  and 
imaginary  output  data  are  required.  By  simply  trans¬ 
lating  the  grating  G  ,  along  the  .\  •  direction  by  one- 
I'ourth  of  a  period  (i.e..  t  2  phase  shift),  the  light  in¬ 
tensity  distribution  along  the  z  axis  becomes 


/(x)  “  J  f(a)(l  +  sin(ii)Ol|da.  lI4) 

which  is,  of  course,  the  sine  transform  of  the  input  f(a). 
Thus  both  the  real  and  imaginary  parts  of  the  complex 
transform  output  can  be  obtained.  We  may  note  that 
the  shifting  of  the  grating  produces  a  ir/2  phase  shift  for 
all  wavelengths.  The  achromatic!  ty  of  the  system  is 
maintained.  This  is  in  contrast  to  other  interferometric 
systems  using  reflective  or  refracting  optics. 

Once  again,  the  intensity  distribution  along  the  z  axis 
can  be  read  out  by  placing  a  thin  beam  splitter  along  the 
X  -z  plane  as  shown  in  Fig.  4.  The  intensity  distribution 
along  the  beam  splitter  is  mapped  onto  a  convenient 
output  plane  outside  the  fringe  volume  using  a  one- 
to-one  telescopic  imaging  system. 

This  simple  implementation  causes  several  problems. 
The  beam  splitter  has  to  be  very  thin;  a  pellicle  beam 
splitter  might,  therefore,  be  used.  Pellicle  beam 
splitters  unfortunately  have  a  tendency  to  vibrate. 
Although  this  is  not  a  fundamental  problem  and  can  be 
solved  using  special  optics,  we  decided  to  use  an  alter¬ 
nate  arrangement  for  the  demonstration. 

The  light  distribution  in  the  output  fringe  volume  is 
equal  to 


/(x.y.x) 


I  +  CO* 


4tF  X 


da. 


If  we  place  a  fourth  grating  with  spatial  frequency  F  at 
2  =  zo  as  shown  in  Fig.  5,  the  grating  will  demodulate  the 
light  back  to  zero  spatial  frequency,  and  the  output  light 
intensity  becomes 


/ix.l) 


1  +  cos 


4tF£o«\ 
f  I 


da. 


(15) 


which  in  effect  gives  the  cosine  transform  of  fia)  sam¬ 
pled  at  spatial  frequency  AwFzo/f.  By  translating  the 


Fir.  4  Mapping  the  transform  output  in  the  fringe  volume  to  a 
convenient  output  plane  using  a  beam  splitter. 


**NSi.*’  :sc 


Fix.  .*>  Sampling  ot  'ipaiial  trequencies  wuh  a  iranslatinx  fourth 
xraimx. 
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Fig.  6.  Decoding  the  transform  output  with  a  tilted  fourth  grating. 


fourth  grating  along  the  z  axis  over  the  entire  fringe  box, 
the  cosine  transform  output  can  be  scanned  out. 
However,  such  a  translation  technique  takes  time,  and 
our  goal  is  real-time  operation. 

We  note  that  Eq.  (15)  is  not  a  function  of  x  and  y. 
The  sampled  cosine  output  can  be  obtained  by  mea¬ 
suring  the  intensity  at  any  point  of  the  output  plane.  If 
we  tilt  the  fourth  grating  toward  the  x  -z  plane  so  that 
it  traverses  the  fringe  volume  as  illustrated  in  Fig.  6, 
different  parts  of  the  grating  along  the  y  direction  will 
be  sampling  the  intensity  at  different  positions  along 
the  z  direction.  The  output  intensity  distribution  be¬ 
comes 


/(y) 


1  +  coa 


'4irayF  tan0] 

I  /  i 


da. 


which  can  be  written  as 

/(«'><•  r  /(a)da  +  C  f(a)  cosluj'alda. 

Jo  Jo 

where 


(16) 


(17) 


ill'  •  (4)rFv  tanio)//. 

and  0  is  the  tilt  angle  of  the  fourth  grating. 

Once  again,  a  biased  output  of  the  cosine  transform 
is  obtained.  However,  the  spectrum  is  now  displayed 
along  the  u)'  direction,  which  is  orthogonal  to  the  a  axis. 
The  cosine  spectrum,  therefore,  lies  along  a  direction 
that  is  perpendicular  to  the  1-D  input. 

If  we  translate  the  fourth  grating  by  one-eighth  of  a 
period  (i.e.,  ±jr/4  phase  shift),  we  obtain 

/lul  l  =  r  fia\da C  /(al  siniui'atrfa 

•/O  Jfi 

or  the  sine  transformation  of  /(a). 

To  perform  complex  Fourier  transformations,  we 
need  both  the  real  and  imaginary  parts.  The  cosine  and 
sine  transforms  can  be  displayed  simultaneously  by 
using  a  split  grating  as  the  fourth  grating  at  the  output 
as  shown  in  Fig.  7.  The  grating  lines  at  the  two  halves 
of  the  grating  are  phase  shifted  by  jr/4.  This  produces 
a  cosine  transform  at  one  side  of  the  output  and  the  sine 
transform  at  the  other. 

We  note  that  the  tilted  grating  decoding  technique 
may  also  be  applied  to  the  single-grating  interferometer 
by  placing  the  tilted  grating  against  the  first  grating 


transversing  the  fringe  volume.  However,  the  presence 
of  spurious  diffraction  orders  may  substantially  degrade 
the  transform  output.  In  addition,  transform  infor¬ 
mation  near  dc  may  not  be  displayed  due  to  bliKkage 
by  the  edge  of  the  tilted  grating. 

When  only  cosine  transforms  are  performed,  the 
input  must  be  either  single  sided  (i.e..  0  S  o  S  /  )  or  an 
even  function  [i.e.,  /(a)  =  /(-a)|  to  avoid  ambiguity. 
With  a  dual-channel  system,  however,  the  input  is  no 
longer  limited  to  being  single  sided  or  even.  The  pro¬ 
cessor  is  actually  performing  the  decomposition 


)  exp)-(uirt)(ia  ■ 


cost  u>a  Ida 


-•  S-J'"' 


sin(u)a)(fa 


.7i 


where  fin)  is  an  arbitrary  real  and  non-negative  func¬ 
tion,  which  in  general  is  a  combination  of  even  and  odd 
parts.  That  is,  the  decomposition  is  not  in  terms  of  the 
cosine  transform  of  the  even  part  and  the  sine  transform 
of  the  odd  part  but  rather  in  terms  of  7r  and  J/,  which 
are  the  real  and  imaginary  parts,  respectively,  of  the 
Fourier  transform.  In  this  operating  mode,  the  input 
can  extend  from  -/  to  /  with  unique  information  on  both 
sides,  and  the  space-bandwidth  product  of  the  proces¬ 
sor  is  in  effect  doubled.  This  doubling  of  the  space- 
bandwidth  product  can  be  considered  to  occur  because 
the  system  acts  as  a  dual-channel  processor  simulta¬ 
neously  performing  the  real  and  imaginary  parts  of  the 
Fourier  transformation.  A  more  detailed  discussion  on 
the  space-bandwidth  product  of  the  grating  interfero¬ 
metric  processor  is  given  in  a  later  section. 


IV.  Experimental  Remits 

An  interferometric  optical  processor  was  set  up  using 
the  configuration  illustrated  in  Fig.  6.  As  the  input,  a 
ground  glass  diffuser  was  backilluminated  by  an  unfil¬ 
tered  xenon  arc  lamp  to  produce  an  incoherent  light 
source,  which  in  turn  was  used  to  illuminate  a  slit  and 
transparency.  To  test  if  the  system  can  produce  the 
correct  Fourier  transform  output,  a  transparency  with 
an  intensity  transmittance  of  1  +  cos(2Tga)  was  used. 
Thus  the  input  is  given  by  fin)  =  rect[(a  -  rto)Ll[l  + 
cos(27r^o)).  The  cosine/sine  transform  outputs  were 
read  out  using  a  1024 -element  linear  detector  array.  In 
Fig.  8(a),  we  show  the  cosine  transform  output  obtained 
with  the  interferometric  processor.  For  comparison, 
we  computed  the  Fourier  transform  of  a  function: 


-C 


Fi|{.  7.  Simultaneous  cosine  and  sine  transformations. 
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Fig.  tj.  Optically  generatad  coaine  and  sine  transforms  vs  the  digitally 
computed  versions. 


/  («)  «  W(n(  rect  'j  [1  +  cos(2)rga  +  0\\  expiifl),  1 19) 

where  rect((«  —  ao)/L]  =  1  for  ao  —  L/2  5  a  <  ao  +  L/2, 
lV(a)  is  an  apodizing  factor  added  to  account  for  the 
nonuniform  illumination  of  the  input,  «q  is  the  offset 
of  the  input  from  the  axis,  g  is  the  spatial  frequency,  d> 
is  the  phase  of  the  input  cosine  function,  and  (f  is  an 
arbitrary  phase. 

The  choice  of  the  grating  position  that  is  designated 
as  cosine  or  sine  transform  is  arbitrary  and  unimportant 
as  long  as  an  exact  w/2  phase  shift  is  introduced  between 
the  sine  and  cosine  transforms  and  a  consistency  is 
maintained.  The  difference  in  the  outputs  caused  by 
different  initial  choices  of  the  grating  position  is  simply 
a  constant  phase.  ( A  similar  situation  arises  in  coherent 
optical  processing  where  a  coherent  reference  beam  is 
introduced  to  measure  the  complex  output  of  the  optical 
processor. )  In  our  experiment,  the  constant  phase  term 
explifl)  was  obtained  by  varying  H  in  the  computation 
until  the  real  part  of  the  zero  order  of  the  digital  Fourier 
transform  matched  that  of  the  optical  output.  In  Fig. 
8(b).  we  show  the  computed  real  part  of  the  Fourier 
transform.  Except  for  a  small  phase  difference,  there 
is  a  very  good  correspondence  between  the  optically 


obtained  and  digitally  computed  results.  The  small 
phase  error  might  have  been  caused  by  irregularity  in 
the  local  frequency  of  the  gratings. 

The  fourth  grating  was  then  translated  by  one-eighth 
of  a  period  to  produce  a  rr/4  phase  shift  on  the  recom¬ 
bined  beams.  The  optically  obtained  sine  transform 
output  is  shown  in  Fig.  8(c).  For  comparison,  the  cor¬ 
responding  computed  imaginary  part  of  the  Fourier 
transform  is  shown  in  Fig.  8(d).  Once  again,  the  optical 
output  matches  the  digitally  computed  result. 

V.  Camparltoii  to  a  Mfcttolson  StoNar  Marfaromotar 

The  operation  of  this  system  can  be  explained  more 
intuitively  by  going  back  to  the  Michelson  stellar  in¬ 
terferometer.  The  Michelson  stellar  interferometer  in 
its  basic  form  is  simply  two  apertures  separated  by  a 
distance  d  as  shown  in  Fig.  9(a).  The  CDSC  is  obtained 
by  measuring  the  complex  visibility  of  the  fringes 
formed  for  different  separations  of  the  apertures.  We 
see  from  Fig.  9(b)  that  we  have  a  similar  situation  with 
the  grating  interferometer.  The  interfering  light  fields 
are  sheared  by  an  amount  Ax  »  z'  tan{siii~UX/)],  where 
F  is  the  spatial  frequency  of  the  grating.  Thus,  similar 
to  the  Michelson  stellar  interferometer,  two  parts  of  a 
wave  front  separated  by  a  transverse  distance  Ax  are 
made  to  interfere,  and  the  resulting  fringe  visibility 
along  the  z  axis,  therefore,  follows  a  Fourier  transform 
relationship  with  the  source  intensity  distribution. 
Unlike  the  Michelson  interferometer  where  the  visi¬ 
bility  for  each  value  of  Ax  must  be  measured  separately, 
the  grating  interferometer  presents  the  interference 
fringes  for  the  different  shear  distances  simultaneously 
at  different  positions  of  z .  This  interferometric  optic^ 


Fix.  9.  Cumparisiin  with  a  Michelson  Stellar  interferometer:  lai 
Michelson  stellar  interferometer;  ibl  simultaneoua  display  of  all 
■.patial  frequent  les  with  the  Kratmg  interferometer:  ici  wavelength- 
dependent  shearing  of  input  wave  front. 
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Fig.  10.  Shifting  of  fringe  volume  for  off-axis  input  point. 


processor  is,  therefore,  basically  an  interferometer  that 
is  capable  of  measuring  in  parallel  the  1-D  CDSC  of  a 
Held  emanated  from  a  homogeneous  incoherent  line 
source  at  far  Held  as  created  by  the  collimating  lens. 

The  spatial  frequency  as  measured  by  a  Michelson 
interferometer  is  given  by  u  =  2xAx/X/,  which  is 
wavelength  dependent.  A  Michelson  interferometer 
can,  therefore,  operate  only  with  quasi-monochromatic 
radiation.  To  make  the  system  achromatic,  a  variable 
amount  of  shear  must  be  introduced  so  that  Ax(X)/X  * 
constant.  The  amount  of  shear  introduced  by  the 
grating  interferometer  varies  with  wavelength,  since  to 
the  first  order.  Ax  *  2F\z,aa  illustrated  in  Fig.  9(c). 
The  spatial  frequency  is  thus  u  =  AirFz/f,  which  is 
wavelength  independent  and  hence  achromatic. 


VI.  Comparlaon  to  Cottoront  Optical  Proeoasing 

Although  the  interferometric  processor  operates  on 
incoherent  inputs,  its  performance  characteristics  are 
very  similar  to  a  coherent  optical  processor.  For  a  co¬ 
herent  system,  the  transform  relationship  can  be  written 
as 


f(pl  "  J'  /(al  exp  da. 

and  for  the  interferometric  optical  processor, 

•  J'  f(a)  cos  (ja 

and 


(20) 


(21) 


f/(z)  “  ^  /(a)  sin  da.  (22) 

The  difference  between  the  two  processors  is  in  the 
scaling  factors  2x/X/  and  AwFIf.  The  scaling  factor  for 
a  coherent  optical  processor  is  inversely  proportional 
to  wavelength,  while  the  scaling  factor  for  the  interfer¬ 
ometric  processor  is  independent  of  wavelength.  Hence 
the  interferometric  processor  is  equivalent  to  an  ach¬ 
romatic  Fourier  transform  system.^-  ’ 

For  systems  that  use  AO  modulators  or  discrete 
sources  (e.g.,  the  vector-matrix  multiplier)  as  input 
devices,  the  space-bandwidth  product  is  typically  input 
device  limited  and  can  be  easily  defined.  For  a  coherent 
optical  processor  and  an  interferometric  optical  pro¬ 
cessor,  on  the  other  hand,  the  space-bandwidth  product 
is  usually  processor-limited  and  is  dependent  on  many 
factors  such  as  the  exact  optical  configuration,  the  //No. 
of  the  transform  or  collimating  lens,  the  aperture  of  the 
lens,  and  the  performance  of  the  lens.  In  the  following. 


we  compute  the  space-bandwidth  product  (SBWP)  of 
the  single-grating  and  the  multiple-grating  interfer¬ 
ometers.  The  approximations  and  assumptions  used 
in  the  computation  are  listed: 

(1)  The  angular  bandwidth  of  the  gratings  in  the 
processor  are  very  wide  (e.g.,  surface  relief  gratings). 

(2)  The  spurious  diffraction  orders  are  not  a  problem, 
and  no  stops  are  used  in  the  system.  (This  assumption 
is  always  valid  for  systems  that  utilize  beam  splitters  to 
demodulate  the  output.  It  is  also  valid  for  the  multi¬ 
ple-grating  system  using  a  tilted  fourth  grating  for 
output  demodulated  if  the  gratings  are  properly  fabri¬ 
cated.  The  effects  of  spurious  diffraction  orders  will 
be  discussed  in  a  later  section.) 

(3)  The  gratings  are  of  sufficient  size  that  the  entire 
field  over  the  full  angular  bandwidth  is  diffracted  by  the 
gratings. 

(4)  The  field  distribution  along  the  optical  z  axis  is 
considered  as  the  output  plane.  [In  actual  system  op¬ 
eration,  it  is  necessary  to  transfer  this  output  plane  to 
a  more  accessible  plane  using  a  beam  splitter  or  a  grat¬ 
ing.  This  analysis  is  directly  applicable  to  a  system  that 
utilizes  a  beam  splitter  to  demc^ulate  the  output  (Fig. 
4).  With  the  use  of  a  tilted  fourth  grating,  however,  the 
performance  characteristic  may  be  degraded  in  the 
decoding  process,  and  the  system  SBWP  may  be  lower 
than  computed.) 

(5)  First-order  approximations  are  used  to  describe 
the  operation  of  the  interferometer. 

For  both  the  single-  and  multiple-grating  interfero¬ 
metric  systems,  except  for  the  change  in  coordinate 
position,  we  have 

/U)  =  ^  f(a)da  +  ^  ^  fia)  cos(uia)(ia, 

where  w  =  AirFzIf. 

Consider  specifically  the  multiple-grating  processor 
as  illustrated  in  Fig.  10.  The  fringe  box  translates  down 
or  up  depending  on  the  position  of  the  input  point  For 
a  point  input  at  a  =  ao,  the  fringe  box  extends  along  the 
z  axis  up  to 


where 


I  ^sin*' 


'  (\/^  -  sin 


tan*'  |  — 


(23) 


corresponding  to  a  cutoff  spatial  frequency  of 

,  IF  [A  2Dao\ 

fT\2  !  I 

For  the  on-axis  point  oq  »  0,  the  cutoff  frequency  is 
a  maximum  and  is  given  by 


(24) 


(-.(O) 


.4  costf 
X/ 


(25) 


where 


«  “  sin*'if .\)  (26) 

is  the  diffraction  angle  -hown  in  Fig.  10.  Notice  that 
the  cutoff  frequency  is  wavelength  dependent  We  also 
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Fig.  11.  Space-bandwidth  products  of  a  single-grating  interfero¬ 
metric  processor,  a  multiple-grating  interferometric  processor,  and 
a  coherent  optical  processor. 


point  out  that,  except  for  the  cosine  factor,  Eq.  (25) 
describes  the  highest  frequency  that  can  be  coherently 
processed  by  a  Fourier  transform  lens  with  aperture 
width  A,  focal  length  /,  and  the  input  placed  a  focal 
length  in  front  of  the  lens.  In  fact,  an  expression 
equivalent  to  that  of  Eq.  (24)  for  such  a  coherent  pro¬ 
cessor  is  given  in  the  Appendix. 

Now  consider  an  input  signal  of  length  21,  extending 
from  a  =  —I  to  a  =  1.  The  end  point  I  establishes  a 
cutoff  frequency  vdD-  Although  one  could  deHne  a 
space-bandwidth  product  (SBWP)  given  by  the  prod¬ 
uct  of  the  line  length  by  the  maximum  cutoff,  i.e., 
4lJ'r(0),  it  is  more  meaningful  for  arbitrary  inputs  to 
establish  the  SBWP  only  as  specified  over  the  uniform 
(nontapering)  part  of  the  spectrum  as  given  by  the 
lowest  cutoff  frequency  i>c(0.  Thus  we  define  the 
SBWP  as 

SBWPd)  =  4/y,in,  (27) 

which  with  substitution  from  Eq.  (24)  becomes 


Equations  (24)-(28)  also  hold  for  the  single-grating 
interferometer  by  letting  D  =  0.  The  single-grating 
interferometer  can  be  considered  to  be  the  limit  of  the 
multiple-grating  interferometer  as  D  approaches  zero. 
It  is  seen  from  Eq.  (28)  that  the  SBWP  is  maximized  for 
D  =  0,  i.e.,  for  the  single-grating  interferometer. 

It  is  worthwhile  to  compare  the  above  results  with 
those  for  a  cla.ssical  Fourier  transform  lens  used  as  a 
coherent  optical  processor  (see  Appendix).  The  SBWP, 
defined  in  a  similar  way  as  for  Eq.  (28),  for  the  coherent 
processor  is 

coh|SBWP(20|  -  2//X  •  (sintun->((/l  -  2/)//|l).  (29) 

where  the  input  extends  from  -  A/2  to  2/  -  A/2. 

In  Fig.  11.  the  SBWP  is  plotted  as  a  function  of  the 
length  of  the  input  for  the  coherent  optical  processor 
and  for  the  two-grating  interferometer  geometries. 
Note  that  for  the  multiple-grating  interferometer  and 
the  coherent  optical  processor,  there  is  an  optimum 


input  length.  To  maximize  further  the  SBWP  of  the 
multiple-grating  interferometer,  it  is  necessary  to 
minimize  the  length  of  the  interferometer  2D.  (This 
is  equivalent  to  moving  the  input  closer  to  the  Fourier 
transform  lens  in  the  case  of  a  coherent  optical  proces¬ 
sor.)  As  we  shorten  the  interferometer,  the  gratings 
eventually  collapse  onto  a  single  plane,  and  we  essen¬ 
tially  have  the  single-grating  interferometer.  (This  is 
similar  to  putting  the  input  right  in  front  of  the  trans¬ 
form  lens  of  the  coherent  optical  processor.) 

The  above  space-bandwidth  product  is  computed 
from  geometric  considerations.  A  second  factor  that 
can  limit  the  space-bandwidth  product  of  the  system 
is  the  dynamic  range  of  the  output  detector  array.  The 
bias  level  is  determined  by  the  integrated  input  signal 
level 

B  =  ^  f(a)<ia. 

As  I  increases  so  does  the  bias  level.  Thus  the  number 
of  input  points  that  can  be  processed  for  the  worst  case 
is  equal  to  the  dynamic  range  of  the  output  detector. 
The  space-bandwidth  product,  however,  can  be  much 
larger  if  the  transform  of  the  input  consists  of  only  a 
small  number  of  points  (e.g.,  a  sinusoidal  input).  In 
that  case,  the  space-bandwidth  product  of  the  system 
is  limited  by  the  system  geometry. 

We  should  note  that  the  same  dynamic  range  limi¬ 
tation  also  applies  to  a  coherent  optical  processor  if  the 
complex  output  is  to  be  read  out  with  a  detector  array. 
However,  there  is  more  flexibility  with  the  coherent 
system.  For  example,  if  the  input  data  are  on  a  carrier 
and  the  information  around  dc  is  unimportant,  the 
reference  beam  can  be  lowered  to  a  level  corresponding 
to  the  highest  signal  level  in  the  desired  portion  of  the 
output  spectrum.  With  the  interferometric  processor, 
on  the  other  hand,  the  bias  is  fixed  at  a  level  corre¬ 
sponding  to  the  dc  output. 

We  emphasize  that  unlike  many  other  optical  systems 
that  utilize  gratings  (e.g.,  the  achromatic  Fourier 
transform  system),  spurious  diffraction  orders  do  not 
necessarily  pose  a  problem  for  the  grating  interfero¬ 
metric  processor.  With  the  use  of  a  beam  splitter  to 
decode  the  output,  the  spurious  orders  have  minimal 
effect,  since  they  propagate  in  directions  away  from  the 
output  plane  as  long  as  the  angular  subtend  of  the  input 
is  less  than  the  diffraction  angle  of  the  grating.  This  is 
true  for  both  the  single-grating  and  multiple-grating 
arrangements  utilizing  beam  splitters.  For  a  multi¬ 
ple-grating  system  that  employs  a  tilted  grating  to  de¬ 
modulate  the  output,  the  spurious  diffraction  orders  are 
also  not  a  serious  problem.  First,  the  higher  diffraction 
orders  propagate  away  from  the  output  plane,  and  for 
the  right  arrangement  they  will  not  pose  problem.  The 
zero  order  does  propagate  toward  the  output  plane,  but 
it  can  be  suppressed  significantly.  For  example,  if  a 
thin  sinusoidal  phase  grating  is  used,  the  powers  in  the 
various  diffraction  orders  follow  the  Bessel  functions. 
By  adjusting  the  exposure  in  the  fabrication  of  the 
grating  to  where  the  zeroth-order  Bessel  function  goes 
to  zero,  the  dc  term  can  be  substantially  suppressed. 
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Furthermore,  the  zeroth  order  is  incoherent  with  the 
first  orders  that  are  being  combined  (large  path  differ¬ 
ences).  Thus  the  effect  of  the  dc  term  that  remains  is 
to  add  a  small  uniform  bias  to  the  output.  Spurious 
orders,  however,  do  create  serious  problems  when  using 
a  tilted  grating  to  decode  the  output  of  a  single  grating 
system,  and  such  an  arrangement  is  generally  not  de¬ 
sirable. 

VII.  AppIleatlOfM 

The  interferometric  optical  processor  can  be  config¬ 
ured  to  perform  operations  currently  accomplished  by 
other  optical  processors.  In  addition,  it  can  perform 
operations  that  cannot  be  easily  done  otherwise,  since 
it  operates  on  incoherent  light. 

For  example,  to  process  a  temporal  signal,  the  scan¬ 
ning  electron  beam  of  a  CRT  can  be  used  to  convert  the 
temporal  electrical  signal  into  a  spatial  formatted  light 
signal.  The  processed  output  is  received  by  a  detector 
array  which  is  read  out  after  integrating  over  the  time 
of  a  complete  scan  of  the  electron  beam.  Such  a  system 
can  operate  over  a  large  range  of  frequencies  by  varying 
the  scan  rate  of  the  electron  beam.  The  ease  in  ob¬ 
taining  the  real  and  imaginary  parts  of  the  complex 
Fourier  transform  output  gives  the  interferometric 
processor  a  potential  advantage  over  optical  systems 
using  AO  m^uktors  and  SAW  devices. 

The  interferometer  is  a  parallel  processing  device,  and 
its  capability  is  not  fully  utilized  when  used  to  process 
single-channel  temporal  signals.  For  applications 
where  parallel  inputs  are  available  (e.g.,  from  an  an¬ 
tenna  array),  the  inputs  can  be  used  to  drive  individual 
LEDs  in  a  linear  array.  Together  with  an  array  of  in¬ 
dividually  addressable  detectors,  the  system  can  be  used 
to  perform  n -point  discrete  Fourier  transforms  at  a  rate 
determined  by  the  frequency  response  limit  of  the 
LEDs.  As  compared  to  the  vector-matrix  multiplier, 
the  interferometric  processor  has  the  advantage  of 
simplicity  requiring  only  simple  gratings  and  collima- 
tion  optics  instead  of  complicated  2-D  masks  and 
beam-forming  optics. 

More  important,  however,  is  the  ability  of  the  inter¬ 
ferometric  optical  processor  to  operate  on  completely 
incoherent  inputs.  The  system,  therefore,  has  the  po¬ 
tential  to  perform  Fourier  transformations  on  signals 
directly  from  the  real  world  using  light  emanating  from 
the  object  scene  itself. 

VIII.  Summary 

We  have  introduced  an  optical  processor  which  is 
capable  of  performing  complex  spatial  1-D  Fourier 
transforms  on  incoherent  optical  inputs.  One  possible 
implementation  of  the  processing  approach  using  a 
grating  interferometer  was  presented,  and  the  operation 
of  an  interferometric  optical  processor  is  experimentally 
demonstrated.  The  potential  ability  of  the  system  to 
directly  process  images  using  the  light  emanated  from 
an  object  scene  removes  limitations  imposed  by  input 
devices  and  opens  new  avenues  for  optical  processing. 
The  research  on  this  type  of  interferometric  processor 
is  still  in  it  initial  stages,  and  much  has  yet  to  be  devel- 
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Fig.  12.  Geometry  of  a  coherent  optical  procesaor  for  the  processing 
of  real  inputs. 


oped.  The  optical  configurations  presented  are  only 
examples  of  many  other  possible  implementations  of  the 
interferometric  optical  processor.  Alternative  system 
configurations  with  different  performance  character¬ 
istics  that  are  optimized  for  specific  applications  can  be 
developed. 


The  work  was  supported  by  the  U.S.  Army  Research 
Office.  Parts  of  this  paper  were  presented  at  the  In¬ 
ternational  Optical  Computing  Conference,  Apr. 
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Appendix:  Space-BandwIdlh  Product  of  a  Coharanl 
Optical  Procaaaor 

For  an  optical  Fourier  transform  lens  to  process  the 
input  signal  without  loss  of  information,  all  the  light 
emanating  from  the  input  aperture  must  be  passed  by 
the  lens  aperture,  i.e.,  no  vignetting.  Consider  the  1-D 
optical  processor  illustrated  in  Fig.  12.  At  position  a 
the  cutoff  frequency  defined  by  geometry  consider¬ 
ations  is  given  by 

«  A"'  sin|taii'‘|(/4  -  2|a|)/2fi|.  (All 

Thus  the  cutoff  frequency  varies  along  the  input  aper¬ 
ture.  It  is  possible  to  define  a  maximum  space-band¬ 
width  product  via 

SBWP„„  -  4/kc(0),  (A2) 

where  the  input  is  assumed  to  extend  from  —Ito+l  and 
the  extra  2  comes  from  the  two  sidebands.  However, 
to  insure  that  the  input  is  processed  without  toss  of  in¬ 
formation,  the  highest  spatial  frequency  of  the  input 
must  be  limited  to  vdl).  Then  the  SBWP  becomes 

SBWPKI  »  4/i'r(n.  (A3I 

If  the  input  is  real,  the  transform  is  Hermitian.  One 
sid'^band  of  the  output  spectrum  can  be  discarded 
without  loss  of  information.  By  placing  the  input 
asymmetrically  from  -A/2  to  21  -  A/2,  the  (one-sided) 
cutoff  frequency  becomes 

I'  (2/1  ”  sinlun'Ml  't  ~  -G/7|l,  iA4l 

where  0  <  /  <  A/2.  The  SBWP  becomes 

SBWP'(/I  «  J/v./J/l.  '.A.tI 
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To  determine  the  optimum  input  size  for  a  given 
value  of  A  and  f,  we  let  21  =■  rxA  and  take  the  partial 
derivative  of  Eq.  (A3)  to  get 
(2n-U 

—  «  —  (A6) 

( 1  -  n  ip 

and  of  Eq.  (A5)  to  get 

(2n-l)  .A-’ 

- -  —  .  (Alt 

1 1  -  n  )■>  p 

For  lenses  with  //Nos.  of  //2  or  slower,  n  is  nearly  0.5  for 
both  cases.  That  is,  the  optimum  SBWP  tends  to  be 
obtained  when  the  input  aperture  size  is  one-half  of  the 
lens  aperture  size. 
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ABSTRACT 

Synthetic  aperture  imaging  with  active  systems  is  relatively  well 
developed.  The  implementation  of  passive  systems  on  the  other  hand, 
is  being  hindered  by  low  SNR  and  the  difficulty  in  achieving  two- 
dimensional  aperture  fill.  The  achromaticity  of  the  grating  inter¬ 
ferometer  permits  the  use  of  a  wider  instantaneous  spectral  band¬ 
width  to  produce  better  SNR  and  it  provides  a  simple  means  to  synthe¬ 
size  a  two-dimensional  aperture.  A  passive  synthetic  aperture  imag¬ 
ing  approach  implemented  with  an  achromatic  grating  interferometer  is 
described  and  laboratory  experimental  results  are  presented  to  demon¬ 
strate  the  proposed  concept. 


I.  INTRODUCTION 


Synthetic  aperture  imaging  is  a  familiar  development  in  the 
microwave  wavelengths  in  the  form  of  active  SAR  for  fine  image  reso¬ 
lution  imaging  without  resorting  to  the  use  of  impossibly  large 
antennae  [1-3].  Relative  motion  between  the  sensor  and  the  object 
scene  is  employed  to  generate  synthetic  aperture  data.  In  parti¬ 
cular,  Walker  showed  [4]  that  through  polar  formatting,  two- 
dimensional  synthetic  apertures  of  rotating  objects  can  be  generated 
with  a  range-Ooppler  imaging  system.  An  imaging  radar  operating  in 
such  a  mode  is  often  referred  to  as  a  spot-light  SAR  since  the  radar 
is  locked  onto  a  fixed  point  in  the  target  scene.  Aleksoff  [5]  has 
also  demonstrated  the  interferometric  imaging  of  rotating  objects  in 
the  visible  wavelength  using  a  fringe  projection  approach.  Synthetic 
aperture  systems  utilizing  either  the  range-Ooppler  or  the  fringe 
projection  approaches  are  active  systems.  That  is,  radiations  are 
launched  by  the  imaging  system  and  the  fields  reflected  by  the  object 
scenes  are  measured.  In  this  paper,  a  passive  approach  implemented 
with  an  achromatic  grating  interferometer  for  fine  resolution  imaging 
of  incoherent  scenes  is  described. 

Passive  synthetic  aperture  imaging  techniques  (or  aperture  syn¬ 
thesis)  have  been  employed  in  astronomical  imaging  for  sometime.  The 
radio  telescope  is  a  prime  example.  An  array  of  antennae  are  used 
with  the  earth  rotation  to  generate  the  synthetic  aperture  [6]. 
While  the  operation  of  the  active  systems  is  based  on  coherent 
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diffraction  theory,  the  operation  of  passive  systems  is  based  on  the 
van  Cittert-Zernike  theorem.  A  coherent  active  system  makes  use  of 
the  Fourier  transform  relationship  between  the  complex  amplitudes  at 
two  widely  separated  (object  and  detection)  planes.  The  synthetic 
aperture  is  generated  by  measuring  the  complex  amplitude  of  the 

reflected  radiation  at  the  detection  plane.  The  passive  approach  on 
the  other  hand,  takes  advantage  of  the  Fourier  transform  relationship 
between  the  intensity  distribution  of  an  incoherent  field  at  the 
object  plane  and  the  mutual  coherence  of  the  field  at  the  detection 
plane.  The  synthetic  aperture  is  therefore  generated  by  measuring 
the  spatial  coherence  of  the  field  at  the  detection  plane.  The 

Fourier  transform  kernel  is  scaled  by  wavelength.  This  fact  permits 
synthetic  aperture  imaging  systems  with  fixed  baselines  to  generate 
radial  fills  of  the  synthetic  aperture.  This  wavelength  dependency 
however,  also  limits  the  instantanteous  spectral  bandwidth  that  can 
be  used  with  either  the  active  or  the  passive  systems.  Too  large  an 
instataneous  bandwidth  can  produce  a  smearing  effect  in  the  phase 

history  signal  (transform  data)  and  reduce  the  system  signal-to-noise 
ratio.  This  restriction  does  not  pose  any  problem  for  active  systems 
but  it  limits  the  performance  of  passive  interferometric  imaging  sys¬ 
tems.  Natural  emissions  and  illuminating  sources  are  spectrally  very 
wide  band  and  the  SNR  of  the  images  reconstructed  with  a  passive 

interferometric  imaging  system  is  directly  proportional  to  the 
instantaneous  spectral  bandwidth.  System  sensitivity  can  therefore 
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be  improved  if  a  means  can  be  found  to  increase  the  instantaneous 
spectral  bandwidth  without  suffering  the  signal  loss  due  to  the  wave¬ 
length  dependency  of  Fourier  transformation. 

In  this  paper,  we  introduce  the  use  of  an  achromatic  grating 
interferometer  for  passive  synthetic  aperture  imaging  of  rotating 
objects  in  the  visible  region.  The  achromaticity  of  the  grating 
interferometer  permits  the  instantaneous  spectral  bandwidth  to  be 
widened  substantially  while  maintaining  the  capability  to  use 
wavelength  diversity  to  achieve  aperture  fill.  And  as  we  shall  show, 
the  achromatic  grating  interferometer  also  offers  other  operational 
advantages  over  conventional  interferometric  systems  when  used  for 
passive  synthetic  aperture  imaging.  A  system  design  approach  and 
preliminary  experimental  results  demonstrating  the  proposed  passive 
synthetic  aperture  imaging  concept  are  presented. 


II.  PASSIVE  INTERFEROMETRIC  IMAGING 


Consider  the  system  geometry  as  illustrated  in  Figure  1.  An 
object  field  A  radiates  according  to  its  blackbody  temperature  and 
the  emitted  radiation  is  received  by  an  antenna  pair  located  at 
and  02*  The  emitted  electromagnetic  field  is  both  temporally  and 
spatially  incoherent.  However  by  filtering  the  received  signals  with 
a  narrowband  filter,  the  field  may  be  treated  as  temporally  coherent. 
(The  filter  bandwidth  necessary  to  achieve  coherency  will  be  dis¬ 
cussed  later.)  Let  £^(3^  6^)  and  £^(“2' 

disturbances  on  the  object  plane.  By  using  the  Huygens-Fresnel  prin¬ 
ciple,  the  wave  disturbances  at  (x^  y^ )  and  {x2,  y2)  due  to 
the  radiation  from  the  object  scene  can  be  written  as  [7] 

-i(l<R^-u)t) 

£(xi ,  y^ )  =  7  //  >  sp  - — ^ ds^  (1 ) 

z  ^ 

and 


-i  {kR2-'jt) 

£(x2,  y2)  *  r  //  ^o^®2’  *^“2  *^^2 

where  \  is  the  wavelength  corresponding  to  the  center  frequency  of 
the  filter,  R^  and  R2  are  the  distance  between  (a^,  3^)  - 
(xp  y^ )  and  (32,  82)  -  (X2.  y2^  respectively,  k  is  the 

wave  number  and  u  the  angular  frequency  of  the  radiation,  and  I 
represents  the  support  of  the  object  scene.  Since  the  radiation  is 
assumed  to  be  quasimonochromatic,  we  shall  from  here  on,  neglect  the 
time  varying  phase  factor  and  let  E  =  E 

0 


The  mutual  intensity  between  the  fields  at  (x^,  )  and 

(xg,  y2)  is  given  by 

*yi )^*( ^2 ’■1^2 ^ ^ 

where  <•>  denotes  ensemble  averaging.  The  complex  degree  of  spatial 
coherence  (CDSC)  is  related  to  the  mutual  intensity  by 


Using  the  Fresnel  approximation,  it  can  be  shown  that  [6] 


SS  I(a,  0)  e  ’2>T(ua+VB) 

V  ^ - 1 - 

//  I(a,  b)  da  dB 
X 

where  u  =  (x^  -  X2)/Rx  and  v  »  (y^  -  y2)/Rx  and 

'i'  =  C  (Xi  +  y^  )  -  (x^  +  y^)]/2RX.  The  CDSC  .  ui2(u,  v  ). 
is  equal  to  the  Fourier  transform  of  the  object  intensity  distribu¬ 
tion  multiplied  by  a  quadratic  phase.  If  the  far  field  condition  is 
satisfied,  then  0  and  e’^’^  =  1,  the  CDSC  becomes 

»12(u,  V)  =  K  //  I(a,  B)  e  (5) 


where  K  is  a  complex  constant. 

We  see  that  the  CDSC  is  related  to  the  intensity  distribution  of 
the  object  field  by  a  Fourier  transform  relationship.  We  can  there¬ 
fore  reconstruct  the  object  image  by  performing  an  inverse  Fourier 
transform  on  the  CDSC  function. 


III.  CONVENTIONAL  APPROACHES 


The  COSC  can  be  obtained  with  Young's  simple  two-pinhole  arrange¬ 
ment  or  equivalently,  Michelson's  stellar  interferometer,  by  measur¬ 
ing  the  complex  visibility  of  the  interference  fringes.  The  complex 
fringe  visibility  function  is  defined  as 

2)  I(x^,y^) 

V(x^,y^;  X2,y2)  =  ~  Kx^y^)  +  irx2.y75 

Except  for  very  near  field  conditions,  we  can  generally  make  the 
approximation,  Ux^y^))  =  I(x2,y2)  and  we  have 

v(u,v)  =  u(u,v).  (3) 

To  measure  the  COSC  via  the  visibility  function  is  not  very  efficient 
since  the  interference  pattern  must  be  scanned  out  with  a  small 
detector  to  determine  the  phase  and  modulation  depth  of  the  sinu- 


soidal 

fringe 

pattern. 

An  alternate 

approach  is  to 

use  a  beam 

splitter 

to 

recombine 

the  fields 

gathered  at 

(x^,y^)  and 

( ^2  ’•1^2  ^ 

such 

that  they 

propagate 

coincidentally. 

The  entire 

output  can  then  be  focussed  on  and  detected  by  a  photo  detector. 
The  output,  as  we  shall  show,  corresponds  to  the  inphase  portion  of 
the  Fourier  transform  sampled  at  spatial  frequency  (u,  v).  The 
quadrature  part  can  be  obtained  by  shifting  the  position  of  the  beam 
splitter  to  generate  a  n/Z  phase  difference  between  the  recombined 
fields.  By  changing  the  positions  of  (x^,y^)  and  i^2'y2'> 

different  spatial  frequencies  of  the  object  is  measured.  A  large 

aperture  can  then  be  synthesized  in  time. 


In  interferometric  imaging  as  described  above,  the  detected  field 
is  assumed  to  be  quasimonchromatic.  Natural  emission  and  illuminat¬ 
ing  sources  tend  to  be  spectrally  very  wide  band.  To  obtain  co¬ 
herency,  the  instantaneous  spectral  bandwidth  of  the  received  signal 
must  be  limited  with  a  narrow  band  spectral  filter.  In  the  follow¬ 
ing,  we  examine  the  maximum  instantaneous  spectral  bandwidth  that 
can  be  used  with  a  conventional  system.  A  two-dimensional  analysis 
will  be  used  for  two  reasons.  First,  the  number  of  terms  is  reduced, 
making  the  analysis  easier  to  understand.  Second,  the  interferometer 
itself  is  a  one-dimensional  system,  performing  at  any  single  moment, 
the  correlation  of  the  fields  between  two  points  in  the  detection 
plane.  A  two-dimensional  model  comprised  of  the  direction  of  pro¬ 
pagation  and  the  direction  of  the  lateral  separation  (or  shears)  of 
the  fields  to  be  correlated,  thus  provides  an  adequate  description 
of  the  system  operation  at  any  one  time. 

Consider  a  single  point  radiator  in  the  far  field  as  illustrated 
in  Figure  2.  The  source  is  assumed  to  be  located  so  far  away  that 
the  field  detected  by  the  two  detection  apertures  can  be  considered 
to  be  parts  of  a  plane  wave  with  the  angle  of  incidence  determined 
by  the  position  of  the  source.  The  fields  at  the  two  apertures  can 
be  described  by 


El(ao)  = 


-ill  [(R+Ar)  ] 


-i~  C{R-Ar)yi-AS‘'] 


=  /o(aQ)  e  (9) 

where  r’  is  the  target  range,  o(a^)  is  the  field  intensity  at  the 
detection  plane  due  to  a  point  source  located  at  a  distance  from 

the  optical  axis,  and  p  =  sin  (  atn  \)  =  Sq/Rx  is  the 

spatial  frequency  of  the  incident  plane  wave.  (For  a  2-dimensional 
object  field,  a^  is  the  position  of  the  point  source  after  it  is 
projected  onto  the  a  axis.) 

The  total  path  difference  from  the  source  to  the  two  detection 
apertures  is 


2ar 

where  S  is  the  separation  of  the  apertures.  I*  the  fields  are  made 
to  interfere  with  each  other,  the  correlation  output  is  equal  to 

I(<Xo)  =  C  o(a^)  [l  -  COs(^  /l  -  xV  )]  (10) 

where  C  is  a  complex  constant.  Using  the  Fresnel  approximation,  we 
have 


[(a^)  =  C  o(a^)  1  -  cos 


(1  - 

2ira  S 

0 

~rs — 

(1  - 

i)11 


For  small  a  , 
0 


we  can  further  approximate  the  correlation  output  intensity  and  obtain 

=  C  oCa^)  [l  +  cos  (aiTUa^)]  (12) 

where 


u 


S 


=  U  • 


Integrating  over  the  entire  object  scene,  we  obtain  a  cosine  trans¬ 
formation  relationship  plus  a  bias  term. 

Iq(‘J)  =  C  /  o(oIq)  da^  +  C  /  o(ajj)  COS  (2^00^)  do^  (13) 


by  adjusting  the  optics  such  that  a  ir/2  phase  difference  is  inserted 
between  the  two  interfering  fields,  it  can  be  shown  that  the  correla¬ 
tion  output  intensity  becomes 

=  C  /  o(a^)  da^  +  C  /  0(0^)  sin  (Z-rua^)  da^  (14) 
Z  Z 

Both  the  inphase  and  quadrature  components  of  the  Fourier  transform 
can  therefore  be  obtained. 


The  transform  kernel  u  =  S/xR  is  wavelength  depen¬ 
dent.  To  see  how  wide  an  instantaneous  spectral  bandwidth  can  be 
used  with  a  conventional  interferometer,  let  us  restrict  the  maximum 
phase  misregistration  to  1/4  of  a  wave.  That  is. 

Sot-  Sa  1 

J  0  — 

xR  ■  (X  ax)R  -  4  • 
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The  instantaneous  spectral  bandwidth  must  then  be  limited  to 


AX  < 


x^R 

-  -t) 


(15) 


for  =  0  (i.e.  an  on-axis  point  source),  we  have  ax  <  x,  the 
spectral  bandwidth  is  unlimited.  This  is  to  be  expected.  For  a 
point  source  at  the  optical  axis,  the  path  difference  to  the  two 
detection  apertures  is  zero. 


The  instantaneous  spectral  bandwidth  decreases  with  increasing 
With  the  "spotlight"  mode  of  synthetic  aperture  imaging,  the 
system  field  of  view  (FOV)  is  defined  by  the  footprint  of  the  real 
aperture  and  it  is  locked  on  to  a  fixed  point  in  the  target  scene. 
Thus,  for  a  point  target  located  at  the  edge  of  the  FOV,  we  have 
a  Rx/A  where  A  is  the  diameter  of  the  detection  aperture.  The 
instantaneous  spectral  bandwidth  that  can  be  used  with  a  conventional 
system  for  passive  synthetic  aperture  imaging  is  therefore  equal  to 


AX  <  — 

4( 


X  R 

3Rx  x^ 
- 5- 


) 


XA 

=  ^  ' 


(16) 


For  example,  if  x  =  0.55  um,  and  A/S  =  1/20,  then  ax  <  6.3  nm. 


IV.  ACHROMATIC  GRATING  INTERFEROMETER 


The  achromatic  grating  interferometer  has  been  studied  exten¬ 
sively  and  applied  to  many  applications  [8-11].  Tai  and  Aleksoff  in 
particular,  has  utilized  it  as  an  incoherent  optical  processor  for 
real  time  Fourier  transformation  [11  J*  The  transform  operation  is 
limited  to  a  single  dimension  and  it  could  create  difficulties  in 
some  applications.  This  limitation  however,  does  not  represent  any 
drawback  for  synthetic  aperture  imaging  since  under  most  operating 
conditions,  a  long  baseline  is  available  only  along  a  single  dimen¬ 
sion  (e.g.,  along  the  fuselage  of  an  airPlsne). 

For  the  same  reasons  stated  earlier,  a  two-dimensional  model  will 
be  used  to  describe  the  operation  of  the  grating  interferometer.  In 
Figure  3,  we  show  the  basic  geometry  of  a  grating  interfermeter  that 
can  be  used  for  passive  synthetic  aperture  imaging.  It  is  composed 
of  three  gratings  of  the  same  spatial  frequency  F.  Light  waves 
impinging  on  gratings  G-j  and  G2  are  diffracted  and  recombined  at 
G^  which  is  placed  a  distance  d  from  the  plane  formed  by  G^  and 
G2.  The  gratings  G^  and  G2  also  represents  the  input  apertures 
and  they  have  a  width  of  A  and  are  separated  by  a  distance  S.  After 
being  diffracted  by  G^,  the  recombined  fields  propagate  coinci¬ 
dentally  and  are  focussed  down  onto  a  photo  detector  of  width  W^. 
Let  us  assume  as  before,  that  the  object  plane  is  sufficiently  far 
away  that  the  light  field  emanated  by  a  point  source  in  the  object 
scene  may  be  considered  to  be  a  plane  wave  as  it  impinges  on  the 
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input  apertures  of  the  grating  interferometer.  Thus,  for  a  single 
point  source  of  intensity  Of  a  )  located  at  a  distance  a  from 
the  optical  axis  as  shown  in  Figure  3,  the  field  right  before  it  is 
diffracted  by  the  grating  G3  is  equal  to 

£(3^)  =  yO(a^)  I  exp  [  i2^(p  -  F)x  -  ikR’  Jl~-  -  ikd  J  1  -  x^(p  -  F)2] 

+  exp  [i2ir(p  +  F)x  -  ikR'  /  1  -~x^p^  -  ikd  /  1  -  x^{p  +  F)^]  | 

(17) 

where  F  is  the  spatial  frequency  of  the  diffraction  gratings,  p  = 
sin  [atn  (■3q/R)]/x  and  k  =  2ir/\.  The  corresponding  intensity  dis¬ 
tribution  can  be  written  as 


Ij,(a  )  =  CO(aQ)  I  1  ^  cos  [4,Fx  -  X'fp  +  F)^  -  [T^ip  -  F7)]j 

(13) 

where  C  is  an  aooropriate  constant.  Using  a  grating  G3  with  spatial 
frequency  F  to  recombine  the  beams  diffracted  by  G1  and  G2,  the  field 
intensity  is  demodulated  back  to  base  band  and  the  intensity  distri¬ 
bution  right  after  being  diffractec  by  G3  becomes 


C0(=.^) 


COS 


X‘(p 


-  £)' 
(19) 
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If  we  expand  the  square-root  terms  into  binomial  series  and  keep  only 
the  first  order  terms,  we  obtain 


ZFda 


^  UU  m 

Ig(aJ  =  CO(a  J  [  1  -  COS  2-(— ^)J 


(20) 


We  see  that  the  output  intensity,  to  the  first  order,  is  not  a  func¬ 
tion  of  wavelength.  The  object  scene  within  the  FOV  of  the  inter¬ 
ferometer  can  be  considered  as  a  superposition  of  incoherent  point 
sources.  The  output  intensity  due  to  the  entire  scene  can  then  be 
expressed  as 


I^(u)  a  C  y'o(a)  da  +  C  y'  0(a)  COS  (2itUo)  do 


(21) 


where 


u 


2Fd 

R 


which  gives  a  cosine  transform  relationship  together  with  a  bias 
term.  If  we  translate  grating  G3  laterally  by  1/3  of  a  grating 
period,  the  phases  of  the  two  fields  being  recombined  are  shifted 
achromatically  by  ±^,'4.  It  can  be  shown  that  the  output  intensity 
becomes 


I^(:j)  3  C  0(a)  da  +  C  /o{^)  sin  (2TUa)  da  (22) 

giving  the  sine  transform  of  the  intensity  distribution  of  the  inco¬ 
herent  source. 

For  a  given  value  of  F,  d  and  R,  the  grating  interferometer  mea¬ 
sures  the  source  spectrum  at  a  spatial  frequency  of  u  =  Fd/R.  Unlike 
the  conventional  system  described  earlier,  the  transform  kernel  is 
not  wavelength  dependent. 
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V.  ACHROMATIC  I TY 


The  achromaticity  of  the  grating  interferometer  can  be  explained 
more  intuitively  by  noting  that  the  spatial  frequency  measured  is 
equal  to  the  u(x)  =  S/xR  where  S  is  the  lateral  separation  of  the 
fields  being  correlated.  To  make  the  measurements  achromatic,  it  is 
necessary  that  S  be  also  wavelength  dependent  such  that  u  =  S(x)/xR  = 
constant.  The  grating  interferometer  satisfies  this  requirement  to 
the  first  order  by  the  fact  that  different  wavelengths  are  diffracted 
at  different  angles  such  that  at  a  given  distance  d  where  the  fields 
from  the  two  inputs  apertures  are  recombined,  the  separation  S  of 
the  fields  being  correlated  are  approximately  proportioned  to  x  as 
illustrated  in  Figure  4.  The  range  of  lateral  shear  for  given  grat¬ 
ing  separation  S  is  limited  to  S  *  A/2.  If  the  interferometer  is 
achromatic,  then 

=  Constant.  (23) 

j  ^  X  / 

Thus, 


R(x  +  ax/2)  R(x  -  ax/2^ 

S+A~S-A 

and  the  instantaneous  spectral  bandwidth  of  an  achromatic  system  is 
equal  to 


AX^=-^  (25) 

The  bandwidth  is  3  times  wider  than  that  of  a  conventional  system  as 
computed  earlier  (see  Eq.  16).  The  grating  interferometer  is 
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achromatic  only  to  the  first  order.  There  is  no  guarantee  yet  that 
it  is  achromatic  over  the  desired  spectral  bandwidth  of  . 

In  the  following,  the  achromaticity  of  the  grating  interferometer  is 
analyzed  using  a  third  order  analysis. 

In  the  previous  section,  we  have  shown  that  for  a  point  source 
located  a  distance  a  from  the  axis,  the  output  intensity  can  be 
approximated  to  the  first  order  by 

lQ(a)  =  C0(a)  ^1  -  cos  ( — ^“)] 

The  grating  interferometer  is  thus  achromatic  to  the  first  order  if 
the  far  field  condition  is  met.  If  we  expand  Eq.  19  and  keep  all 
the  terms  up  the  third  order  however,  we  obtain 

1^(3)  2  CO(oi)  I  1  -  [cos  2TXpFd  (2  +  +  x^p^  +  |  F^x^  +  |  p^F^x'* 

(27) 

Since  the  FOV  is  confined  to  within  the  foot  print  of  the  real  aper¬ 
ture,  p  generally  is  small  and  we  can  further  aoproximate  the  output 
intensity  distribution  by 

I^(a)  =  C0(a)  {  1  -  cos[^^(  2  *  x^f2  +  ^4^)]}  (28) 

We  find  that  the  system  is  not  achromatic  to  the  third  order  even  if 
the  object  plane  is  located  at  far  field.  While  the  field  is  demod¬ 
ulated  achromatically  (i.a.  the  spatial  frequency  of  the  output  light 
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intensity  distribution  is  zero),  the  phase  of  the  output  cosine  func¬ 
tion  is  wavelength  dependent.  Let  us  once  again,  limit  the  maximum 
misregistration  to  with  1/4  wave.  That  is. 


Fdo 

IT 


+ 


(x 


ax)2  f2  . 


Then,  the  spectral  bandwidth  of  the  system  must  be  limited  to 


ax2=  R/4aFd(2xF^  +  |  f\^)  (30) 

For  an  on-axis  point  source  (a  =  0),  the  system  is  completely  achro¬ 
matic.  The  system  however,  must  remain  achromatic  for  all  points 
within  the  FOV.  For  a  synthetic  aperture  imaging  system,  the  FOV  is 
defined  by  the  footprint  of  the  real  aperture.  Thus,  the  worst  case 
represents  a  =  Rx/A.  And  since  d  =  aS(x)/2  tan  [9(x)]  where  aS(x) 
is  the  shear  distance  and  ©(.x)  is  the  diffraction  angle  for  the 
center  wavelength  x,  the  spectral  bandwidth  can  be  rewritten  as 

ax2=  Ax/2as(x)  cos  e(x)j  2[sin^  8(a)]/\  +  SCsin"^  9(A)]/2x|  (31) 

We  see  that  the  achromatic  spectral  bandwidth  is  highly  dependent  on 
the  diffraction  angle  e  (or  equivalently,  the  soatial  frequency  of 
the  gratings).  To  increase  the  achromatic  bandwidth,  ©  must  be  made 
small.  However,  reducing  the  diffraction  angle  will  increase  the 
overall  size  of  the  interferometer.  The  light  waves  from  the  two 
input  apertures  have  to  travel  a  longer  distance  before  they  overlap. 
Thus,  the  optimum  spatial  frequency  for  the  gratings  in  the 
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interferometer  is  one  that  provides  an  achromatic  bandwidth  that  is 
just  sufficient  for  the  system  operation. 


Since  the  range  of  shear  distances  is  limited  by  the  system  geo¬ 
metry  to  AX^=2AX/S,  the  grating  interferometer  only  has  to  be  achro¬ 
matic  over  AX^.  The  optimum  spatial  frequency  for  the  gratings  in 
the  interferometer  can  be  obtained  by  making  the  third  order  achro¬ 
matic  bandwidth  in  Eq.  31  equal  to  the  desired  achromatic  bandwidth 
in  Eq.  25.  The  two  spectral  bandwidths  are  equal  when 


4  cos  9(x) 
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or  e(x)  a  20.5’.  For  a  center  wavelength  of  x  =  0.55  um,  the  optimum 
spatial  frequency  for  the  gratings  will  be 


12L±(±i  »  700  Ip/mm. 

If  a  FOV  larger  than  2Rx/A  is  desired  (a  larger  FOV  can  be 
achieved  by  simply  using  a  larger  detector  since  the  FOV  is  defined 
by  P  *  where  P  is  the  point  spread  function  of  the  real  aperture 
projected  out  to  the  object  plane,  is  the  width  of  the  photo 
detector  and  *  denotes  convolution,  then  a  lower  spatial  frequency  F 
must  be  used  to  maintain  the  achromatic  property  over  the  entire  FOV. 


VI.  APERTURE  FILL 


For  a  given  interferometer  geometry,  the  correlation  output  pro¬ 
vides  the  value  of  the  object  spectrum  at  a  single  spatial  frequency. 
The  system  geometry  must  be  changed  to  generate  transform  data  at 
other  spatial  frequencies  in  order  to  synthesize  a  two-dimensional 
aperture.  First,  if  there  is  relative  rotational  motion  about  the 
optical  axis  between  the  object  scene  and  the  interferometer,  trans¬ 
form  data  can  be  gathered  along  a  circular  arc  in  the  transform 
plane.  Second,  if  either  the  baseline  or  the  detection  wavelength 
is  varied,  transform  data  along  a  radial  line  in  the  transform  plane 
is  obtained.  Combining  these  two  effects,  a  two-dimensional  aperture 
can  be  synthesized.  The  means  to  achieve  such  an  aperture  fill  will 
be  discussed  below. 

Let  us  assume  for  now  that  the  interferometer  is  looking  straight 
down  at  the  object  scene  located  at  the  optical  axis.  The  interfer¬ 
ometer  output  provides  the  transform  data  at  a  single  point  in  the 
Fourier  transform  plane.  As  the  object  (or  equivalently,  the  inter¬ 
ferometer)  rotates  about  the  optical  axis,  data  is  gathered  along  a 
circular  arc.  To  see  what  the  phase  history  will  look  like  as  the 
target  rotates  relative  to  the  sensor,  consider  a  single  point  target 
located  at  a  distance  from  the  optical  axis.  Its  Fourier  spec¬ 
trum  is  a  linear  phase  in  the  x-direction  and  the  in-phase  or  quad¬ 
rature  portion  of  the  Fourier  spectrum  is  a  sinusoidal  pattern  with 
spatial  frequency  of  approximately  a^/Rx.  The  interferometer 
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measures  the  Fourier  spectrum  of  the  object  scene  at  spatial  fre¬ 
quency  u  =  s/Rx  =  2Fd/R,  V  »  0.  As  the  target  or  the  interferometer 
rotates,  the  cosine/sine  transforms  of  the  object  spectrum  is  at  u  = 
cos(e)  s/Rx  and  v  =  sin(0)  s/Rx  are  measured  where  o  is  the  angle  of 
rotation.  The  correlation  output  of  the  interferometric  sensor  as  a 
function  of  time  (phase  history)  can  be  written  as 

(4it  sin  (tan~^  (S/2R))  a  [1  -  cos  (ut  -$  )] 

- r - 

(3i) 

which  can  be  rewritten  as 

lQ(aQ,  $  ,  t)  =  1  +  cos  {  ka^Cl  -  cos  (wt  )]|  (34) 

where  $  is  the  angular  position  of  the  point  target  relative  to  the 
a  axis,  u  is  the  angular  rotation  rate  and  k  is  a  constant. 

In  Figure  5,  we  show  the  correlation  output  as  a  function  of  time 
for  point  targets  located  at  different  positions  (a^,v  )  in  the 
object  scene  as  the  object  is  rotated  relative  to  the  interferometer 
at  a  constant  rate.  Since  an  incoherent  input  scene  can  be  con¬ 
sidered  as  an  superposition  of  incoherent  point  radiators,  the  cor¬ 
relation  output  for  an  extended  source  as  the  object  rotates  can  be 
expressed  as 

N 

^0^^^  ®  Z  {  1  ^  cos  [ka^  (l  -  COS  (-.Jt 

n=l 
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Radial  fill  of  the  synthetic  aperture  can  be  achieved  with  the 
use  of  baseline  diversity  or  wavelength  diversity.  Since  the  spatial 
frequency  measured  is  determined  by  S/Rx,  it  can  be  changed  by  vary¬ 
ing  S,  the  separation  of  the  detection  aperture  (baseline)  or  x,  the 
detection  wavelength.  Baseline  diversity  can  be  accomplished  with 
the  grating  interferometer  by  moving  the  three  gratings  closer  or 
apart  as  shown  in  Figure  6.  However,  it  may  not  the  desired  mode  of 
operation  for  synthetic  aperture  imaging  since  it  requires  the  use 
of  large  input  windows  to  accomodate  the  change  in  baseline.  It  is 
often  more  preferrable  to  keep  the  input  aperture  fixed  and  achieve 
the  radial  aperture  via  the  use  of  wavelength  diversity.  On  the 
surface,  it  may  seem  that  the  use  of  wavelength  diversity  with  an 
achromatic  system  is  contradictory.  We  note  however,  that  the  inter¬ 
ferometer  has  a  finite  instantaneous  bandwidth  and  the  whole  band 
can  be  slid  back  and  forth  to  obtain  the  wavelength  diversity  for 
radial  aoerture  fill.  This  is  accomplished  by  changing  d,  the  sepa¬ 
ration  between  gratings  G1-G2  and  G3.  Recall  from  Eq.  21  that  the 
spatial  frequency  measured  by  the  interferometer  is  equal  to  2Fd/R. 
Changing  d  thus  varies  the  spatial  frequency  measured.  It  can  be 
seen  from  Figure  7  that  changing  d  is  also  equivalent  to  changing 
the  center  wavelength  of  measurement  x.  The  instantaneous  achromatic 
bandwidth  remains  2Ax/S  for  any  d.  The  factor  of  8  gain  in  SNR  is 
maintained  over  the  entire  synthetic  aperture. 


For  a  single  point  source  located  at  (oq,  ^  ),  varying  d  at  a 
constant  rate  will  generate  a  correlation  output  of 


(36) 


where  T  is  the  linear  translation  rate  of  grating  G3  and  k'  is  a 
constant.  The  correlation  output  as  a  function  of  time  is  a  sinusoid 
whose  frequency  is  dependent  on  the  position  of  the  point  source 


The  two  means  of  aperture  (circular  and  radial)  fill  can  be  used 
together  to  generate  a  two-dimensional  aperture.  Unlike  tomographic 
imaging,  synthetic  aperture  imaging  systems  generally  gather  trans¬ 
form  data  over  an  angular  range  much  less  than  90  degrees  and  the 
two-dimensional  aperture  synthesized  is  typically  a  segment  of  the 
an  annulus  as  shown  in  Figure  8.  The  image  (half-power)  resolution 
in  the  ground  plane  is  approximately  equal  to 
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where  ae  is  the  angle  of  rotation.  The  resolution  in  the  a  direction 
I  is  proportional  to  x"  instead  of  x  because  the  size  of  the 

synthetic  perture  is  proportional  to  ax/x.  The  aperture  is  larger 
for  smaller  x  if  ax  is  fixed.  For  example,  a  real  aperture  imaging 
I  system  operating  at  0.5  urn  with  a  10  cm  aperture  can  achieve  a  ground 

plane  resolution  of  5  m  x  5  m  at  a  range  of  1000  Km.  A  passive  syn¬ 
thetic  aperture  imaging  system  with  the  same  real  aperture  size  and 

I 

I 
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an  aperture  separation  of  2  m  on  the  other  hand,  can  achieve  a  ground 
plane  resolution  of  0.63  m  x  0.58  m  at  the  same  distance  with 
=  0.5  urn,  ax  =  0.2  urn,  and  ae  =  25*. 

To  synthesize  a  two-dimensional  aperture  using  the  grating  inter¬ 
ferometric  sensor,  there  are  two  possible  approaches  that  can  be 
taken.  First,  the  grating  G3  can  be  dithered  back  and  forth  along 
the  optical  axis  (changing  d)  as  the  target  rotates.  This  creates  a 
two-dimensional  synthetic  aperture  that  is  continuous  in  the  radial 
direction  as  shown  in  Figure  9(a).  Alternately,  multiple  demodula¬ 
tion  gratings  (G^)  placed  at  different  positions  along  the  z-axis  of 
the  interferometer  can  be  used  to  measure  several  spatial  frequencies 

simultaneously.  As  the  object  scene  rotates,  an  aperture  composed 
of  continuous  arcs  is  obtained  as  illustrated  in  Figure  9(b).  The 
advantage  offered  by  the  former  approach  is  that  only  a  single  demod¬ 
ulation  grating  and  detector  is  required.  However,  in  some  operating 
modes  where  the  time  available  for  data  acquisition  is  restricted,  it 
may  not  be  feasible  to  dither  the  the  grating  fast  enough  to  synthe¬ 
size  an  acceptable  aperture.  The  latter  method  gather  multiple 
transform  data  simultaneously.  Its  drawback  is  the  need  for  multiple 
sets  of  hardware  such  as  gratings,  detectors,  amolifiers  etc.  The 
approach  of  choice  will  be  dependent  on  the  operating  requirements  of 
the  particular  application. 


VII.  IMPULSE  RESPONSE 


In  the  above  formation  for  image  resolution,  the  synthesized 
aperture  was  for  simplicity  assumed  to  be  completely  filled  and  has 
rectangular  dimensions  of  an.,  A^.  However,  the  synthesized 

aperture  is  likely  to  be  a  segment  of  an  annulus  and  it  is  contin¬ 
uously  filled  along  either  radial  lines  or  circular  arcs.  For  a 
single  point  target,  the  shape  of  the  impulse  response  (IPR)  may  not 
be  important;  the  main  concern  is  the  width  of  the  IPR  which  defines 
image  resolution.  In  the  imaging  of  an  extended  object  however,  the 
shape  of  the  IPR  (sidelobe  structure)  plays  a  vital  role  in  the  image 
quality.  In  this  section,  we  examine  what  the  impulse  responses  look 
like  with  different  aperture  fills.  As  described  in  the  last  sec¬ 
tion,  a  two-dimensional  aperture  can  be  synthesized  by:  (1)  dither 
ing  the  grating  G3  along  the  optical  axis  of  :he  interferometer  or, 
(2)  using  multiple  demodulation  gratings  placed  at  different  posi¬ 
tions  along  the  optical  axis.  The  first  approach  produces  continuous 
radial  fill  and  the  later  approach  produces  continuous  circular  arcs. 

First,  we  show  in  Figure  10,  the  transform  of  a  single  radial 
line  obtained  by  translating  linearly  the  grating  G3  along  the  opti¬ 
cal  axis  while  the  object  is  stationary.  Without  object  or  inter¬ 
ferometer  rotation,  high  image  resolution  is  obtained  only  along  the 
a  direction  through  the  use  of  wavelength  diversity.  If  there  is 
object  rotation,  then  a  two-dimensional  synthetic  aperture  can  be 
generated  by  dit  ering  the  grating  G3  back  and  forth  along  the 


_ 

optical  axis  of  the  interferometer.  The  aperture  fill  obtained  by 
dittering  the  grating  is  shown  in  Figure  11  a  and  the  resulting 

impulse  response  is  given  in  Figure  lib.  (The  synthetic  aperture 

used  to  obtain  Figure  11b  contained  25  radials  lines.  Much  fewer 

lines  are  shown  in  Figure  11a  in  order  to  illustrate  the  structure 

more  clearly. ) 

If  the  position  of  the  grating  is  fixed,  then  a  circular  arc  of 
transform  data  as  shown  in  Figure  12a  is  obtained  by  rotating  the 
object  or  the  sensor.  The  bow  tie-like  impulse  response  of  such  an 
aperture  is  shown  in  Figure  12b.  For  this  case,  high  image  resolu¬ 
tion  is  obtained  in  the  a  direction  where  a  large  aperture  is  synthe¬ 
sized  through  the  rotating  motion  of  the  object  or  the  sensor.  By 
using  multiple  sets  of  gratings  and  detectors,  a  two-dimensional 
aperture  fill  as  shown  in  Figure  13a  is  obtained.  The  corresponding 
impulse  response  of  a  synthetic  aperture  containing  25  arcs  is  given 
in  Figure  13b.  The  impulse  responses  obtained  with  the  two 

approaches  for  aperture  fill  are  different  in  fine  structures  but 

their  over  all  image  qualities  are  about  equivalent.  We  may  note 

also  that  the  two  approaches  can  be  combined.  A  small  set  of 

grating/detectors  may  be  used  to  obtain  simultaneous  measurements  of 
several  spatial  frequencies  and  the  gratings  may  also  be  dithered 
together  to  fill  in  the  gaps.  Combining  the  two  approaches  would 
reduce  the  number  of  grating/detector  sets  required  and  decrease  the 
amount  of  grating  travel  needed  to  fill  the  desired  aperture  when 


VIII.  EXPERIMENTS 


A  grating  interferometer  was  set  up  in  the  laboratory  and  exper¬ 
iments  were  performed  to  demonstrate  some  of  the  concepts  described 
above.  It  is  difficult  to  perform  optical  synthetic  aperture  imaging 
experiments  in  the  laboratory  because  of  the  problems  posed  by  the 
short  target  range  that  is  available.  First  of  all,  the  grating 
interferometric  sensor  requires  the  far  field  condition  to  be  satis¬ 
fied  between  the  object  and  detection  planes.  Moreover,  image  reso¬ 
lution  at  optical  wavelengths  at  such  short  ranges  are  very  high. 
This  would  require  the  use  of  either  a  very  small  sensor  or  a  very 
small  target  scene. 

The  gratings  we  used  in  the  experiments  were  fabricated  on 
dichromated  gelatin  using  two  plane  waves  in  a  symmetric  geometry  to 
produce  a  spatial  frequency  of  500  lo/mm.  Iseally,  grating  G1  and 
G2  should  be  faoricated  with  an  asymmetric  geometry  such  that  dif¬ 
fraction  favors  a  single  order  while  the  grating  G3  should  be  fabri¬ 
cated  with  a  symmetric  recording  geometry  to  produce  two  (±1)  orders 
of  equal  intensity.  However,  to  ensure  that  511  the  gratings  have 
identical  spatial  frequency,  the  same  recording  geometry  was  used 
for  al  1  the  gr ati ngs . 

The  experiment  was  performed  on  a  4'  x  8'  optical  bench.  It  pro¬ 
vided  a  target  range  of  about  1.5  m.  In  order  for  the  FOV  to  be  of 
workable  size,  a  real  collection  aperture  of  0.5  mm  was  used.  The 
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size  of  the  corresponding  point  spread  function  at  the  object  olane 
was  1.5  mm.  The  apertures  were  separated  by  10  mm,  giving  a  A/S 
ratio  of  1:20.  Since  the  aperture  separation  is  so  small,  gratings 
G1  and  G2  were  implemented  by  masking  off  a  single  grating  instead 
of  using  two  individual  gratings.  To  simulate  the  far  field  condi¬ 
tion  at  such  a  short  target  range,  an  achromatic  lens  with  a  focal 
length  of  1.5  m  i''  placed  in  front  of  the  grating  as  shown  in  Figure 
14  to  collimate  the  incoming  field.  To  study  the  responses  of  the 
grating  interferometer  and  compare  it  with  theory,  a  rotating  wide¬ 
band  point  source  was  generated  by  rotating  a  pinhole  that  is  back 
illuminated  by  an  unfiltered  Xenon  arc  lamp.  To  be  a  "point"  source 
to  the  sensor,  the  size  of  the  source  must  be  smaller  than  the  period 
of  the  spatial  frequency  being  measured.  For  the  system  geometry 
used  in  the  experiment,  it  means  that  the  source  diameter  must  be 
less  tiian  75  ym.  A  25  ym  pin  hole  was  used  in  our  experiments. 
Grating  G3  was  placed  at  an  nominal  distance  of  17.5  mm  from  G1-G2 
and  it  was  mounted  on  a  stepping  motor  driven  linear  translation 
stage  in  the  Z  direction  and  a  piezoelectric  driven  translator  in 
the  X  direction.  The  output  light  field  was  focussed  on  a  fiber 
optics  probe  and  detected  by  a  photomultiplier.  The  outout  intensity 
was  then  digitized  and  transferred  to  a  computer. 

In  Figure  15a,  we  show  the  sensor  output  with  the  geometry  of  the 
grating  interferometer  fixed  and  the  point  source  rotating  over  an 
angular  range  of  *45°.  it  is  identical  to  the  phase  history 


described  by  Eq.  34  and  presented  in  Figure  5.  The  grating  G3  was 
then  translated  by  1/3  of  a  fringe  producing  a  ±tt/4  phase  shifts  on 
the  recombined  fields.  The  quadrature  data  was  obtained  as  shown  in 
Figure  15b. 

To  observe  the  ability  of  the  system  to  resolve  images  in  the  8 
direction  with  the  use  of  object  rotation  to  build  a  synthetic  aper¬ 
ture,  the  point  source  was  placed  at  a  starting  position  as  shown  in 
Figure  16a  and  the  output  is  compressed  to  form  a  point  image  in  the 
8  direction.  In  Figure  16b,  we  show  the  ideal  theoretical  sensor 
output  and  in  Figure  16c,  the  experimentally  obtained  data.  The 
corresponding  compressed  image  is  shown  in  Figure  16d.  Next,  the 
starting  point  of  the  rotator  point  source  was  retarded  to  simulate 
a  second  point  target  off-set  along  a  circular  arc  as  shown  in 
Figure  16e.  The  sensor  output  is  given  in  Figure  16g  and  the  com¬ 
pressed  image  in  Figure  16h.  We  see  that  the  compressed  image  is 
off-set  accordingly. 

To  observe  image  compression  in  the  x  direction  using  wavelength 
diversity  to  generate  the  synthetic  aperture,  the  grating  G3  was 
translated  at  an  uniform  rate  while  the  point  target  was  kept  sta¬ 
tionary.  In  Figure  17a,  we  show  the  position  of  the  point  source 
and  in  Figure  17b,  the  sensor  output.  It  is  a  sinusoidal  pattern  as 
described  in  Eq.  36.  The  compressed  image  is  shown  in  Figure  17c. 
The  point  source  is  then  moved  closer  to  the  optical  center  and  the 
resulting  sensor  output  is  shown  in  Figure  17d  and  the  corresponding 


compressed  image  is  given  in  Figure  17e.  Once  again,  the  difference 
in  the  location  of  the  two  point  sources  along  the  a  direction  is 
resolved. 

The  experiments  performed  were  simple,  demonstrating  separately 
aperture  synthesis  along  the  two  dimensions  using  object  rotation  and 
wavelength  diversity.  Nevertheless,  the  experimental  results  demon¬ 
strate  very  well  the  feasibility  of  using  an  achromatic  grating 
interferometer  for  passive  synthetic  aperture  imaging  of  spectrally 
wide  band  object  scenes.  Work  is  continuing  toward  the  demonstration 
of  passive  synthetic  aperture  imaging  of  two-dimensional  extended 
objects. 


IX.  COMPARISON  WITH  CONVENTIONAL  INTERFEROMETERS 


We  have  shown  in  Section  V  that  a  subtantially  wider  instan¬ 
taneous  spectral  bandwidth  can  be  used  with  the  grating  interfer¬ 
ometer  than  with  a  conventional  interferometer,  thereby  giving  the 
grating  interferometer  better  sensitivity.  There  are  also  opera¬ 
tional  advantages  in  using  the  grating  interferometer.  In  a  conven¬ 
tional  system,  the  measurement  wavelengths  and  thus  the  spatial  fre¬ 
quencies,  are  selected  with  a  bank  of  narrow  band  spectral  filters. 
With  the  grating  interferometer  on  the  other  hand,  it  is  determined 
by  the  position  of  the  gratings  which  can  be  varied  continuously. 
Therefore,  continuous  radial  fill  is  difficult  to  achieve  with  con¬ 
ventional  interferometers  (without  changing  the  baseline)  but  simple 
to  implement  with  a  grating  interferometer.  Due  to  the  sampling 
requirements,  measurement  at  many  wavelengths  is  often  required  to 
generate  an  adequate  aoerture  fill.  The  number  of  parallel  channels 
required  could  be  fairly  large  (>30).  The  grating  interferometer 
provides  the  option  of  using  a  single  or  a  small  number  of  parallel 
channels.  The  aperture  fill  is  obtained  by  simoly  dithering  a  single 
or  a  small  set  of  gratings. 

The  grating  interferometer  is  also  less  sensitive  to  system 
vibrationsand  accurate  phase  shifts  can  be  more  easily  achieved 
since  the  phases  of  a  diffracted  waves  are  determined  mainly  by 
the  lateral  position  of  the  grating  fringes.  Moving  a  grating  later¬ 
ally  by  one  fringe  period  (5  urn)  produces  one  full  wave  of  phase 


X.  SUMMARY 


In  this  paper,  we  have  introduced  the  use  of  a  grating  interfer¬ 
ometer  for  passive  synthetic  aperture  imaging.  Since  natural  emis¬ 
sions  and  illuminating  sources  are  spectrally  very  wide  band,  the 
larger  instantaneous  spectral  bandwidth  permitted  by  the  grating 
interferometer  provides  the  grating  based  sensor  with  better  system 
sensitivity  than  conventional  systems.  Moreover,  the  grating  inter¬ 
ferometer  provides  a  simple  means  to  achieve  radial  aperture  fill 
through  the  use  of  wavelength  diversity  and  the  interferometer  is 
more  tolerant  to  system  vibration.  Possible  operating  modes  that 
can  be  used  to  generate  two-dimensional  synthetic  apertures  were  pre¬ 
sented  and  the  resulting  impulse  responses  were  compared.  Simple 
experiments  were  performed  to  demonstrate  the  concepts  proposed  and 
the  results  were  consistent  with  theories. 
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Figure  11:  Impulse  Reponse  of  a  Synthetic  Aperture  Generated  by  Dittering 

Demodulation  Grating  G^.  a)  Synthetic  Aperture  and  b)  Reconstructed 
Image  of  Point  Target.  g1 


Experimental  Results  Demonstratlny  Image  Compression  in  the  |3-Direction 
Using  Object  Rotation  for  Aperture  Synthesis. 
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gure  17:  Experimental  Results  Demonstrating  Image  Compression  in  the 
a-Direction  using  Wavelength  Diversity  for  Aperture  Synthes 


